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Abstract 

Linear perturbations of Friedmann-Robertson- Walker universes with any curvature and cosmo- 
logical constant are studied in a general gauge without decomposition into harmonics. Desirable 
gauges are selected as those which embody best Mach's principle: in these gauges local inertial 
frames can be determined instantaneously via the perturbed Einstein field equations from the dis- 
tributions of energy and momentum in the universe. The inertial frames are identified by their 
'accelerations and rotations' with respect to the cosmological frames associated with the 'Machian 
gauges'. In closed spherical universes, integral gauge conditions are imposed to eliminate motions 
generated by the conformal Killing vectors. The meaning of Traschen's integral constraint vectors 
is thus elucidated. For all three types of FRW universes the Machian gauges admit much less 
residual freedom than the synchronous or generalized harmonic gauge. Mach's principle is best 
exhibited in the Machian gauges in closed spherical universes. Independent of any Machian mo- 
tivation the general perturbation equations and discussion of gauges are useful for cosmological 
perturbation theory. 
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I. INTRODUCTION 



Einstein preferred a finite universe, bounded in space, over an infinite one because he 
wanted to avoid posing boundary conditions. What Einstein really disliked was that in 
open universes some of the motion of inertial frames is due to dragging by matter while 
the rest is due to the boundary conditions at infinity. Thirty-four years ago the authors 
of acclaimed "Gravitation" [1] commented on the Einstein view in a footnote on p. 704: 
"Many workers in cosmology are skeptical of Einstein's boundary condition of closure of 
the universe, and will remain so until astronomical observations confirm it" . The Wilkinson 
Microwave Anisotropy Probe has now provided data [2] which, among many other things, 
have constrained the present value of total mass-energy density parameter of the universe 
to be f2 = 1-02 ± 0.02. With such a result, all three basic sets of standard Friedmann- 
Robertson- Walker (FRW) cosmological models (see, e.g., [1], [3]) are compatible: the models 
with flat spatial sections (with curvature index k = 0,Q Q = 1), positive spatial curvature 
models (k = +I,f2 > 1), as well as negative curvature (k = — l,f2 < 1)- Nevertheless, the 
WMAP data 'marginally prefer' k = +1 (see, in particular, [4]), and, indeed, recently several 
authors studied closed models again in detail (see, e.g., [5], [6]), after years of preference of 
flat universes which have been considered as natural outcomes of inflation. Even an idea 
going back to the Eddington-Lemaitre cosmology has now been revived: if our universe is 
closed today, it was always closed, and perhaps inflation is 'past-eternal' - the universe, 
dominated at early times by a single scalar field, could have started asymptotically from 
an initial Einstein static universe which enters an inflationary expanding phase, succeeded 
then by standard evolution (see [7] and references therein). The recent growing evidence for 
the existence of a cosmological constant A has been an inspiration for the reconsideration 
of spatially closed universe of de Sitter type [6]. 

In the present work we do not, technically, bestow a privilege to any value of spatial 
curvature. All three cases k — 0, ±1 are analyzed in equal detail, and we even discuss, albeit 
briefly, closed hyperbolic and closed flat universes with multiconnected topologies. From the 
physical (to some extent perhaps 'philosophical') point of view, we adhere to the Einstein 
preference, i.e., to the closed universes with standard (spherical) topology, because our work 
on cosmological perturbation theory has been motivated by Mach's principle. 
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A. Mach's principle 



Mach's principle has acquired certain unpopularity among some relativity and cosmology 
circles. The primary reason is perhaps the fact that under that name a range of meanings and 
interpretations, sometimes even mutually contradictory, has gradually accumulated. During 
the Prague conference in 1988 to celebrate the 150th anniversary from Mach's birth [8] 
and, in particular, at the Tubingen conference in 1993 devoted entirely to Mach's principle, 
numerous interpretations have been given (see the excellent book [9]). More recently, Bondi 
and Samuel [10] listed the 'zeroth' plus ten other versions of Mach's principle and described 
within which theoretical framework a particular statement of the principle applies - see also 
[11], where the main formulations from [9] and [10] are summarized. A brief history of Mach's 
principle and its meaning in general relativity and cosmology is given in the Introduction to 
our first paper on the subject [12]. 

Despite a possible scepticism as regard the role of Mach's principle in contemporary 
cosmology, most of the standard treatises on the subject do include a discussion of the 
principle (see, e.g., [1], [3], [13], [14]), and no one can deny Mach's ideas have been a source 
of inspiration to many, not only Einstein. One of the purposes of the present work is to 
demonstrate that a search for a framework in which Mach's principle can be best embodied 
in the cosmological perturbation theory can lead to practical results, such as the formulation 
of the perturbation theory in a completely general gauge, followed by the selection of an 
advantageous 'Machian gauge' for solving specific problems. 

What then do we mean by Mach's principle? As in our previous work [12], as a starting 
point we adapt the original Bondi's formulation from his classical 'Cosmology' [15]: "Local 
inertial frames are determined through the distributions of energy and momentum in the 
Universe by some weighted averages of the apparent motions". More specifically, we turn 
primarily to those among Einstein's equations for linear perturbations of the FRW models 
which represent the constraints, i.e., under suitable conditions partial differential equations 
of elliptic type, connecting the 'initial values' for matter perturbations with the perturbations 
of the metric. In [12] we went quite a way in realizing the Machian program. We studied the 
frame-dragging effects due to slowly, rigidly rotating, but collapsing or expanding spheres 
in the (inhomogeneous) Lemaitre-Tolman-Bondi universes, and we analyzed the dragging 
effects of vector perturbations of the FRW universes described in a special gauge such that 
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three (momentum) constraint equations enabled us to determine instantaneously metric 
perturbations h° k (k = 1,2,3) in terms of energy-momentum perturbations 5T k . In the 
open universes, these are determined uniquely by requiring the perturbations to vanish at 
infinity - rotations are 'absolute' in this sense. In closed universes a linear combination of six 
Killing vectors (three rotations plus three quasitranslations) may be added to the h° k . We still 
obtain the solutions of the three constraint equations when angular momenta corresponding 
to the three rotations and quasimomenta corresponding to the three quasitranslations of 
the sources (determined by 5T k ) are given. In this sense no absolute rotations exist in 
closed universe, only differences of rotation rates are determinable - in accord with Mach's 
ideas that 'all motions are relative'. If, however, the velocities of the bodies, described 
by perturbations of perfect fluid, are given, the metric perturbations h° k are determined 
uniquely. 

The last result is related to the fundamental fact that six globally conserved quantities, 
corresponding to the six Killing vectors in a FRW universe, must all vanish if considered for 
the whole closed universe. The conserved quantities, being the derivatives of superpotentials, 
can be expressed as surface integrals like an electric charge by using Gauss' theorem. As 
the volume surrounded by the surface is expanded over all the universe at a given time, 
the surface must shrink to zero. It was, among others, an attempt to understand Mach's 
principle in cosmological perturbation theory, which inspired us to formulate conservation 
laws with respect to curved backgrounds [16]. The resulting 'KBL superpotential', using the 
designation by Julia and Silva [17], was found, after applying certain natural criteria, to be 
unambigous and most satisfactory in spacetimes with or without a cosmological constant, 
in any spacetime dimension D > 3 (see [17], [18]). It also found applications in the recent 
studies of the causal generation of cosmological perturbations seeding large-scale structure 
formation and of the back reactions in slow-roll inflation (see [19], [20], and references 
therein). 

In the present paper we study general linear perturbations of the FRW universes from a 
'Machian perspective'. This leads us to investigate both rotations and accelerations of local 
inertial frames in perturbed universes, and to develop all the perturbed Einstein equations 
in a general gauge 'ab initio'. 
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B. Cosmological perturbation theory 

Observational evidence for isotropy and homogeneity of the Universe shows that it is 
broadly well described by a FRW model, but the clustering, the galaxies and the stars 
constitute local perturbations from the idealized substratums of cosmological models. The 
other goal of cosmological perturbation theory is to link the physical conditions in the 
early universe with structures observed today From the pioneering work of Lifshitz in 
1946 (see review in [21]), there appeared numerous papers studying linear perturbations of 
FRW universes; for the more recent extensive reviews, see, for example, [22], [23], [24], [25], 
[26]. Recently, several authors have even found the impetus and skill to write down the 
complicated system of equations for the second-order perturbations of the FRW models - 
see [27] and references therein - but applications of these are yet to appear. 

In any cosmological perturbation theory, two problems at once confront us: (i) What 
smooth cosmological model is best suited to our Universe, and (ii) how do we map the 
points of our inhomogeneous Universe onto a chosen cosmological model. Both these prob- 
lems are, in relativity, connected with the gauge freedom that changes the apparent form 
of the perturbations by which our Universe differs from the smooth cosmological model 
adopted. Although the first problem is primarily connected with a difficult question of an 
appropriate averaging of an inhomogeneous universe, or so-called 'fitting problem' (see, e.g., 
[28]), it is also related to the gauge problem because two different FRW universes can be close 
to one another, at least for some time, so one can be considered as a perturbation of another. 
When we have chosen an appropriate background, then we are confronted with the freedom 
to choose coordinates in the real Universe differently and so to remap the Universe onto the 
background model. This gives rise to the commonly discussed gauge freedom. Of course, 
what happens in the real world is independent of what background is used and how we map 
onto it. This is the motive behind gauge invariant perturbation theory, and this made the 
work of Bardeen [29] in 1980, in which gauge invariant quantities combined from cosmolog- 
ical perturbations were first introduced, so influential. The gauge problem is explained in 
technical terms from different (equivalent) perspectives in depth in the literature: one can 
consider a one-parameter family of 4-dimensional manifolds, with Mq a background and M £ 
a perturbed universe, embedded in a 5-dimensional manifold N and connected by a "point 
identification map" which is specified by a vector field X on N transverse to the manifolds 
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M; the gauge transformation is then a change of X (see, e.g., [23], [30]). Alternatively, in a 
more physical vein, in any chosen coordinate system in the real Universe one assigns to all 
physical quantities Q(x a ) also their background values Q(x a ). These, in contrast to Q(x a ), 
do not change their functional dependence on coordinates under an infinitesimal coordinate 
transformation (see [24]). Mathematically, any of these approaches lead to the changes of 
physical quantities as they appear in the following (see, in particular, Section IV). 

There exists a well-known lemma [30] stating that the linear perturbation of a quantity 
is gauge invariant only if the quantity vanishes on the background or is a constant (scalar 
or linear combination of products of Kronecker delta). The density perturbation 5p, for 
example, is not gauge invariant since p is time-dependent function in the FRW backgrounds. 
That is why to obtain a gauge invariant quantity, one has to consider, e.g., the gradients of 
density perturbation or combine 5p with some other quantities. However, as in the black- 
hole perturbation theory, solving for gauge invariant quantities may not mean finding all 
quantities of interest. For example, in the problem of the motion of a charged black hole 
in a weak asymptotically uniform electric field, there is only one gauge invariant quantity. 
We need to fix the gauge at the end in order to find all perturbations of the metric and 
electromagnetic field to see how the hole accelerates [31]. It is advantageous at the start to 
have the possibility of a gauge choice according to the problem in hand. Selecting a gauge 
which implies a physically preferable coordinate system may eventually give both a better 
physical understanding and an easier mathematical procedure. After all, motions in the 
solar system can be described as seen from a frame that rotates to keep a planet 'at rest', 
but are much more readily comprehended in Newton's inertial axes. 

Last but not least, physical effects associated with Mach's ideas like the dragging of iner- 
tial frames are of a global nature and they do require the introduction of suitable coordinate 
frames (the 'gauges'). A true understanding of inertia and inertial frames must involve 
specific frames or coordinates. To borrow Dieter Brill's comment from Ref. [9]: "Mach's 
principle may point the way toward giving physical meaning to quantities usually considered 
frame-dependent" .We return to this issue in the concluding remarks where some of our other 
work on Machian effects [32], [33] is summarized in the context of the present paper. 
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C. The goal 



There exist many frameworks for treating linearized perturbations of FRW universes. 
The one which has been used most frequently involves the synchronous gauge, with all 
quantities decomposed into suitable harmonics in accordance with the spatial curvature. In 
what follows we make a general study of advantageous gauges without imposing a priori 
conditions on them and without decomposition into harmonics. We identify desirable gauges 
as those which embody best Mach's principle. We find that these gauges are also motivated 
by the gauge choices used in full nonlinear general relativity. Most importantly, however, 
they are distinguished by the simplifications they bring both to the perturbed Einstein 
equations and to their physical interpretation. 

What do we mean by Mach's principle within this broader framework? We again start 
from Bondi's formulation that "local inertial frames are determined through the distribu- 
tions of energy and momentum in the Universe by some weighted averages". However, to 
determine a local inertial frame in a general situation means to find both its 'rotation and 
acceleration' from the distributions of energy and momentum, represented by 5Tff. In a 
general situation we need to know the full spacetime metric in a neighborhood of a point 
in order to determine completely local inertial frames at that point; see, e.g. [1], [3]. In 
Wheeler's conception of Mach's principle (e.g. [1], [34]), we have to specify the initial data on 
a Cauchy spacelike surface like the conformal three-geometry and the mass-energy currents, 
solve for the spacetime geometry g^, and thus determine local inertial frames. In general, 
gravitational waves will globally also contribute to the dragging of the inertial frames but 
only when the waves are nonlinear perturbations of a FWR universe. However, limiting our- 
selves to the linear perturbations of the FRW universes, it is interesting to see what data are 
needed to determine the 'accelerations' and 'rotations' of local inertial frames with respect 
to what we call the cosmological (observers') frames. In a perturbed FRW universe, can a 
gauge be found such that the distribution of 5Tj^ determines uniquely and instantaneously 
the rotations and accelerations of local inertial frames via Einstein's field equations? 
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D. The outline 



After first reviewing the properties of a general congruence of timelike worldlines in a 
general spacetime (see, e.g., [35]), we consider the congruence of 'cosmological observers' in 
a perturbed FRW universe with coordinates {x^} as a 'perturbation' of the congruence of 
fundamental observers in the FRW background. We assume that in a 'cosmological gauge' 
{x 11 } the cosmological observers move along x l = constants, but we describe the properties 
of their congruence by covariant expressions which can be calculated in any coordinates. 
The cosmological observers are equipped with their local frame vectors; the timelike ones 
are their 4-velocities, and the spacelike ones lie along their connecting vectors. Now a 
cosmological observer is, in general, accelerated with respect to a local freely falling inertial 
frame, in particular, the one which at a given spacetime point moves with the same 4- 
velocity. Expressing this acceleration in terms of the metric perturbations, we find that only 
certain components of the metric perturbations are needed. Next, we determine the rotation 
of the axes of the cosmological observer with respect to the nonrotating rigid orthogonal 
axes (gyroscopes held in their centres of mass) of the local inertial frame. Having the 
acceleration a. and the angular velocity u> expressed, we have determined the local inertial 
frame: it accelerates and rotates with respect to the corresponding cosmological frame with 
just the opposite vectors, —a and —uj. All these issues are analyzed in Section II. 

Assuming a general congruence of cosmological observers, i.e., equivalently, a general 
gauge, we find that in order to determine the accelerations and (avaraged) rotations of local 
inertial frames in the sense just described, we need to know the metric perturbations <5g o 
and 5g 0i and their first derivatives. The main issue in Sections III and IV is to find and study 
the gauges in which these quantities can be determined instantaneously from the knowledge 
of energy-momentum distributions 5Tjf. We give the perturbed Einstein equations for all 
three types of FRW universes with any value of A, in an arbitrary gauge. We first adopt 
the 'relativists' attitude' and start from the perturbed FRW metric in the form 

ds 2 = (g^ + h^) dx»dx v = dt 2 - a 2 (t)f kl dx k dx l + h^dx^dx", (1.1) 

where the spatial background metric is fki(x m ),k,l,m... = 1,2,3; t is the 'cosmic time', so 
Sg^ u = h^ u . Perturbations h^ v are small so that quadratic terms can be neglected. In one 
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of the standard coordinate systems the background FRW metric g reads 



ds 2 = dt 2 - a 2 {t) 



dr 2 



+ r 2 (d6 2 + sin 2 6d^ 2 ) , 



(1.2) 



1 — kr 2 



where in a positive curvature (closed) universe (k = +1) r G< 0, 1 >, in flat (k = 0) and 
negative curvature (k = — 1) open universes r G< 0, oo), G< 0,7r >,</? G< 0,27r). We 
shall also employ other common alternatives as, e.g., hyperspherical coordinates, 



with Sfc = sinx, x,sinh \ for, respectively, k = +1,0,-1. The perturbations 5T^ are left 
general, but a perturbed perfect fluid is considered as an example. In Appendix A we give 
all the perturbed Einstein's equations and the Bianchi identities starting off from (1.1); in 
Section III we give them using conformal time 77 and metric perturbations defined as is usual 
in the cosmological literature, e.g., in [22], [24], [26] - again in a completely general gauge. 

We do not decompose the perturbations in harmonics nor do we first separate them 
into the scalar, vector and tensor parts (used, e.g., in [23]). Although both methods are 
very useful in cosmology, they involve nonlocal operations. In order to make Fourier-type 
analyses in the space variables, one needs to know quantities in the whole space, which is 
not 'typical' in cosmology. The splitting of a local perturbation into some scalar, vector 
and tensor perturbations is also nonlocal. Imagine a trivial (zero) perturbation in a given 
domain O, and extend it to an annulus A so that it is nonvanishing there. Hence, in O 
the trivial perturbation will split into nontrivial (scalar, vector, tensor) pieces which depend 
on the extension into A. Therefore, a perturbation which is the sum of scalar, vector, and 
tensor parts cannot be uniquely expressed in terms of the Bardeen gauge invariant variables 
[29] which are defined separately for each part. Without using harmonics or splittings, the 
perturbed Einstein field equations are in a form suitable for searching for solutions in terms of 
Green's functions. How the Green's function approach can reveal new aspects of cosmological 
perturbation theory has been recently indicated by Bashinsky and Bertschinger [36]. 

In Section IV, the main purpose is to motivate and describe geometrically several gauges 
in which the accelerations and rotations of the local inertial frames follow instantaneously 
from the field equations. We call these Machian gauges. We also clarify the residual gauge 
freedom that these gauges admit, and make a comparison with two typically non-Machian 
gauges - the synchronous gauge and the generalized Lorenz-de Donder ('harmonic') gauge. 



ds 2 = dt 2 - a 2 [d X 2 + T? k {d9 2 + sin 2 6d V 2 )] , 
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The Machian gauges turn out to admit much less residual freedom. The freedom represented 
by the gauge transformations generated by the conformal Killing vectors in closed (spherical) 
universes is removed by the integral gauge conditions which we impose. In closed hyperbolic 
universes our Machian local gauge conditions fix coordinates uniquely. 

Finally, in Section V we give the field equations in the Mach 1 gauge and show how 
they can be solved to give the local inertial frames when the distribution of the matter 
energy-momentum is given. We also discuss Traschen's integral-constraint vectors [37] , [38] 
restricting possible STff. According to Traschen and others [37], [39], their existence has 
implications for the Sachs- Wolfe effect and for microwave background anisotropies. Traschen 
considered these vectors in the synchronous gauge. By contrast, in the Mach 1 gauge, these 
constraints become a straightforward consequence of the constraint equations and acquire 
a simple, lucid meaning. We find integral constraints also on quantities not considered 
by Traschen. In closed universes these integral constraints are satisfied automatically as 
a consequence of our integral gauge conditions by which motions generated by the spatial 
conformal Killing vectors are eliminated. In Section V we also list all Green's functions 
known in literature which solve the constraint equations needed for the determination of 
the local inertial frames; some are still unknown. We then review our recent work [40], [41] 
on vorticity perturbations of FRW universes and study their effect on local inertial frames. 
As a second example, we consider perturbations of potential type for which the vorticity 
vanishes. At the end we analyze the 'Machian' question of how uniquely local inertial frames 
are determined in perturbed universes. 

In Concluding Remarks (Sec. VI) we briefly summarize the results and discuss global 
aspects of Mach's principle. In Appendix C the Killing and conformal Killing vectors in 
FRW universes are listed, and those harmonics which are needed in the main text are given. 
In Appendix D we discuss briefly the field equations in the other gauges considered in 
Section IV. 
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II. THE ACCELERATION AND ROTATION OF LOCAL INERTIAL FRAMES 



A. The congruence of cosmological observers 

Consider a general spacetime with coordinates {x^ 1 } in which a congruence of timelike, 
non- intersecting worldlines of "cosmological observers" is given by 

x» = x^y\p) i = 1,2,3. (2.1) 

The choice of fixed y % determines the worldline of a particular observer; p is a parameter 
along the worldline, commonly chosen as either the cosmological time x° = t or the observer's 
proper time r. The cosmological observers use their 4-velocity as their normalized timelike 
frame vector, 

r 11 = (dx^/dr) yi = t»/{g a(3 t a t^, = (dx»/dp) yi . (2.2) 

For spatial frame vectors a cosmological observer naturally takes three independent vectors 
specified by Sy l pointing from him to three other observers of the congruence, orthogonal to 

6af ± = P»5x v = (5» - T»T v ){dx»/dy% 5y\ (2.3) 

As a triad of spatial vectors e^, any three linearly independent vectors proportional to 

5x^_ can be taken. A triad based on the connecting vectors is given at a fixed spacetime 

00 

point and can be extended along the observer's worldline because connecting vectors are 
Lie-propagated (see, e.g., ([42], [43]) along the congruence. This gives 

P>^y = t*5x» ± . (2.4) 

Three independent connecting vectors define the triad of unit spacelike vectors m^y. 

Sx» = 5l {i) m^, m (l)fl m^ = -1, (2.5) 

with no summation over index %. Equation (2.4) implies the propagation equations for scalar 
distances 51^ - the 'generalized Hubble's law' (admitting a possibly anisotropic expansion) 
[35] - and the propagation equations for triad m^y Decomposing the derivative of a 4- 
velocity in the standard manner (e.g. [1], [35]), 

r^-u = r v a^ + uj^ + + \6P^ V , (2.6) 
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the acceleration a^, vorticity (antisymmetric), shear a^ u (symmetric) and expansion 6 
are given respectively by 

a„ = t^t\ (2.7) 

^ = ^P u \t k;X -t X;k ), (2.8) 

<V = \P K ,P^ K .,x + r A;K ) - ±6P, V) (2.9) 

= t£. (2.10) 

We obtain, successively, 

Pfal-X = (w\ + ^ + \eP^)5x\. (2.11) 

^(%) / % = (<V + ±MV)™S) TO M' (2.12) 
^H);,^ = + < + [cj^m^m^] mfo. (2.13) 

B. Cosmological observers in a perturbed FRW universe: The frames 

Consider first an unperturbed FRW model described by metric (1.1) with = 0. 
Fundamental (cosmological) observers move along the worldlines x l = constants with 4- 
velocity r M = (1,0,0,0). As the spatial triad, they take three independent vectors 
perpendicular to r M . These need not be chosen to be necessarily mutually orthogonal if, for 
example, coordinates are used in which f ki in (1.1) is not diagonal [as, e.g., in (C3)]. In 
standard coordinates in FWR backgrounds like in Eqs. (1.2) and (1.3), fki is diagonal and 
the vectors 

^ = (0,0, e (0 „ = (0,fcj (2.14) 
are orthogonal. It is easy to normalize them: 

™& = {-9u)-l [0,C1 , = [0,9iJ , (2-15) 

again with no summation over i, index m — 1, 2, 3. The quantities (2.7)-(2.10) characterizing 
the congruence of the fundamental observers are well known: = ZJ^ U = a^ u = 0, 6 = 3a/a, 
the dot is d/dt. 

In a linearly perturbed FRW universe the metric is given by Eq. (1.1); the indices of the 
first-order quantities (including h^ v ) are shifted by the background metrics g , respectively 
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^ . The congruence of cosmological observers will, in general coordinates, be given by (2.1). 
The frame vectors can be written in the form = r M + 5t^, = + Se^ similarly for 
covariant components, and for o+, u^ u , cr^ u , and 9. In general coordinates these quantities 
can be found easily from the expressions given in Section II A. In the following we shall 
assume that coordinates {x^} represent a "cosmological gauge", in which the congruence of 
cosmological observers is given by x % = y l = constants. We find r M to be given by 

^ = ^ + ^ = (1-1/100,0,0,0). (2.16) 

The spatial triad, determined by connecting vectors orthogonal to r M and lying along coor- 
dinate lines, is 

from which the corresponding unit spacelike vectors can be found: 

™(0 = m (i) + H) = [~ h ^ W - \hu/g u )} • (2.18) 

We gave here both the background and perturbed frames for completeness. In the fol- 
lowing we shall often use just the background frames because these are only needed when a 
small, first-order quantity is projected. 



C. The acceleration of local inertial frames 



We shall designate the local frame of a cosmological observer (CO) given by tetrad r M , 
e^, respectively m^, by the COF - cosmological observer frame. This frame moving along 
x % = constants is, in general, accelerated with respect to local freely falling inertial frames. 
Among the inertial frames there is a frame which, moving at a given spacetime point with 
4- velocity r M , is momentarily at rest with respect to the COF; such a frame is called the 
LIF - local inertial frame 1 . The 4-acceleration of the COF with respect to the LIF is given 
by Eq. (2.7). Using r M given by Eq. (2.16) and perturbed metric (1.1), we find 

a x = (0,a l ), (2.19) 

1 There are of course infinitely many LIFs moving with the 4-velocity r M at a given point. However, they 
differ just by purely spatial transformations or constant shifts of time. Among them, there is also such a 
LIF that its origin coincides with that of a corresponding COF and its acceleration is —a 1 . 
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where 

a 1 = 5T l 00 = g lm (-±h 00 , m + h 0m ). (2.20) 

We see that only hoo and ho m are needed in determining the acceleration of the COF with 
respect to the LIF or, equivalently, the acceleration of the LIF with respect to the COF 
(which is —a 1 ). Spatial metric perturbations do not even enter in the frame components of 
the acceleration because the unperturbed spatial triad is needed to the zeroth order only: 

a(») = e ( i )x a x = e(i }X a x (2.21) 

(similarly with projections on unit vectors m^). Although we calculated the acceleration 
in coordinates adapted to cosmological observers, it is given by a covariant expression (2.7) 
which can be expressed in any coordinates. The result is also invariant under gauge trans- 
formations (see Section IV) since in the background this acceleration vanishes. 



D. The rotation of local inertial frames 



Next we wish to determine the rotation of the axes of the COF with respect to the non- 
rotating rigid orthogonal axes (gyroscopes held in their centres of mass) of the LIF at a 
given point and thus, vice versa, the rotation of the LIF with respect to the COF. 

First consider a cosmological observer carrying a gyroscope described by a spacelike vector 
W^ 1 , perpendicular to r M . The gyroscope is transported along observer's worldline by Fermi- 
Walker transport. Another gyroscope, carried by an inertial observer moving with the same 
r M at a given point, does not rotate relative to W^. However, a vector S^, perpendicular to 
r M , which is transported along the worldline of CO in a general manner, will rotate relative 
to by (D f S^)At, where D F S^ is the Fermi- Walker time derivative defined by (see Fig. 1, 
and e.g. [44]) 

D 9^ 

= P%S? v t v = S%r v + (a„S"V, (2.22) 

where a u is the acceleration (2.7) and S^t^ = was used. For the gyroscope, D F W^/dr = 0. 

Now regarding Eq. (2.11) we see that the left-hand side (l.h.s.) is just equal to the 
Fermi- Walker derivative of the connecting vector so that 

= K„ + < + \6PZ)5x\. (2.23) 
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t + dr 




(D F S")dT 




T 







FIG. 1: The Fermi- Walker time derivative DpS^ (based in part on Fig. 2 in [44]). The cosmological 
observer o with four-velocity carries with himself a gyroscope, represented by the spatial vector 
(dashed arrow), and a spatial vector which are both perpendicular to and identical at 
observer's proper time r. After dr, both and remain perpendicular to but the generally 
propagating vector 5 M will differ from the nonrotating, Fermi- Walker transported gyroscope by the 
Fermi- Walker time derivative (DpS^)dr. 

Therefore, since the congruence of cosmological observers has, in general, a non-vanishing 
vorticity and shear, the connecting vectors rotate with respect to gyroscopes. The last term 
in Eq. (2.23) is proportional to 5x^_an.d represents only a dilation of the connecting vector 
due to the (isotropic) expansion of the congruence. Similarly, unit vectors of the COFs 
rotate with respect to gyroscopes according to Eq. (2.13): 



Turning now to the perturbed FRW universes we find, using from Eq. (2.16) and the 
perturbed metric (1.1), the vorticity (2.8) to have a simple form 




(2.24) 



dr 



LUkl — SuJki — - 



{hok,i — h i,k), u 0a — 5uo a — 0, 



(2.25) 



the shear (2.9) turns out to be 

&ki — b^ki = -hki 
and the expansion (2.10) reads 




(2.26) 




TO 



'm 




a 



(2.27) 
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Since both and cr£ are of the first order in h^ v , on the right-hand side (r.h.s.)of Eq. 
(2.24) only enters and the equation takes the form 



Clearly, the vector rotates relative to the gyroscopes and, hence, a gyroscope will 
rotate relative to COF not only due to a nonvanishing vorticity but also due to the presence 



discussion in [45]); it is not true, as sometimes stated ([26], p. 334) that a spin (a gyroscope) 
precesses relative to the cosmological frame at a rate given just by the vorticity uo k i- 

The axes of a LIF are determined by three orthogonal gyroscopes, while those of a COF 
are determined by three approximately orthogonal vectors or, after their normalization, 
by unit vectors m^. As a consequence, on average the rotation of the COF relative to 
the LIF (moving with the same r M at a given point) is determined just by vorticity u k i. 
Indeed, there is a significant difference between the terms on the r.h.s. of Eq. (2.28): oo k i is 
antisymmetric while a k i is symmetric and traceless. If at a given instant a vector lies 
along a principial direction of o~ k i, its direction will be changed only by the vorticity. As in 
fluid kinematics [46], it is just the vorticity which describes the "effective angular velocity" 
of the fluid (see Fig. 2). 

Hence, we conclude that in order to determine the averaged rotations of local inertial 
frames with respect to the cosmological frames in the perturbed FRW universes, it is sufficient 
to determine the vorticity tensor (2.25), i.e., spatial gradients of h Qk - 

Rather than by uj k i the rotation is usually represented by the "cosmologist's vorticity 
vector" 



where m' 




(2.28) 



of a shear. A gyroscope will precess in a gravitational wave described by h k i (cf. the 



A _ _~\o\xv— 

2 ' 



(2.29) 



where 




(2.30) 
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FIG. 2: Because of the shear with the principal axes a, b, almost orthogonal unit vectors mm, 
m(2) change their directions but they do not, on average, rotate with respect to the fixed axes x 1 , 
x 2 . 

However, curl?; represents twice the effective rigid local angular velocity of the fluid in an 
inertial frame [46]. Therefore, the averaged rigid angular velocity of COFs with respect to 
LIFs is determined exactly by u l given by Eq. (2.31). Equivalently, LIFs rotate with respect 
to COFs with angular velocity — u l . 

E. A note on generalized backgrounds 

When the background is not FRW but, say, a Lemaitre-Tolman-Bondi or Bianchi model, 
the accelerations and the averaged rotations of the LIFs with respect to the COFs can still 
be determined from just h 00 or h i0 . For example, if the coordinates can be chosen such that 
9/j.u = '§ (lv + h^ u , where g 0i = 0, g 00 is an arbitrary function of time x°, and g ik are arbitrary 
functions of x x , then cosmo logical observers given by x % = constant have their accelerations 
with respect to LIFs equal to a 1 = g lm (— |/ioo,m + ^om,o _ !<7 0<> 9oo,o^mo)) an d their vorticity 
is uj^ v = 0, LOki = 5uki = !(<?oo)~ 5 (^ofc,« ~ hoi t k), u 0a = 0. 
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F. Sources and their description in the cosmological frame 



In the FRW universes the background energy-momentum tensor is commonly taken to be 
the perfect fluid stress tensor, = (p + p)U fl U 1 ' — p5 v ^ so that in the comoving coordinates 

Tj = p, T] = -p5], T° = T^ = 0. (2.32) 

The energy density p and the pressure p of the matter can describe a standard perfect fluid 
with a given equation of state. Alternatively, one may regard these expressions as the stress 
tensor components of a homogeneous time-dependent scalar field $ of an inflationary model 
with the energy density p = p$ and effective pressure p = p<$> (see, e.g., [47]). The special 
case with p + p = 0, p = —p = constant, corresponds to the de Sitter vacuum spacetime 
with a cosmological constant A = — p = ~p, commonly interpreted as a vacuum energy. 
Any of these background matter contents can be considered in the present work. We shall 
thus not, in general, specify the form of the perturbations 5T^ of the energy-momentum 
tensor. Employing the frame vectors given by (2.17), we find the frame components of 
perturbations (indicated by [e] and [m])for a general energy-momentum tensor to be given 
by 

^(0) = 5T M(0) = *C 5T tk = (-0«) iW M(O = 5T i ~(P + (2-33) 

^[il(o) = (-<y-^ ]( o) = M$, Km = (9u/9 kk )^T$ ]{l) = ST?, (2.34) 

with no summation over i, k. By employing the "mixed" tensorial coordinate components 
of perturbations, we see that their values, except for 5T®, coincide - up to the background 
"normalization" factors ~ (— ^ ii ) ± i - with their frame scalar components. 
In the case of a perfect fluid the coordinate components read 

5T° = 5p, 5T? = (p + p)(h t0 + Vj, 5T* = (p + p)V\ 5T t k = -5p5f, (2.35) 

where 5p and 5p are perturbations of the matter density and pressure. The velocity 

W = £, (2.36) 

is the spatial part of the perturbation of the fluid's 4- velocity 

U" = U* 1 + 5U» = (1 - I /ioo, V*) . (2.37) 
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It is easy to see that the 4- velocity is approximately the unit timelike vector since we assume 
V 4 C 1, and terms proportional to V 2 and Vh can thus be neglected. The 4-acceleration 
of the fluid is defined by N 1 = U.^,U U . Since the background value = 0, N 1 is of the first 
order. The standard condition A^U^ = thus implies A^U ^ = 0, and hence A = 0. The 
calculation of the spatial components yields 

A* = v* + 2-V 1 - \g is h m , s + g is h s0 = V l + 2-V 1 + a\ (2.38) 

where a % is the acceleration (2.20) of the cosmological frame with respect to the local inertial 
frame, or —a 1 is the acceleration of the LIF with respect to the COF. The acceleration (2.38) 
is the fluid's acceleration with respect to the LIF, whereas V 1 characterizes its acceleration 
with respect to the COF. If the fluid is momentarily at rest in the COF, V % = 0, and the fluid 
has the same acceleration with respect to the COF as the LIF has, V 1 = —a 1 , then A 1 = 0, 
as it should. Since the fluid's acceleration vanishes for the background, its frame components 
are just A^ = (— g u )~^A^ = A\ It is not difficult to check that the acceleration (2.38) 
satisfies the perturbed relativistic Euler's equations, 

(p + p)A i = <5(rr%J , IT" = g ifl - U l U^ (2.39) 

where LT M is the projection tensor into the 3-space orthogonal to U^. As we shall notice in 
Section III B, these are just the spatial parts of the perturbed Bianchi identities. 
The vorticity of the fluid is defined by 

n « = le a ^%Q lS , ^ = 1 -{U T ^ S - U 5 .^). (2.40) 

Since £l a U a = WVL^s = and fi's are of the first order, we again get VL° = = il 0a . The 
nonvanishing spatial parts turn out to be 

tiki = \ l(V k + ZiofcXi ~(Vi + h ol ), k ] , SV = \e m {Vi + h ol ), k = \e ik % k + J- (2.41) 

w 1 is the vorticity vector (2.31) of the cosmological frame. Again, this result is plausible in 

the following sense: Since the LIF rotates with respect to the COF with — uj\ then if the 

fluid rotates with respect to the COF with \e lkl Vi : k = —uj\ it does not rotate with respect 

to the LIF, = 0. As with the acceleration, the frame components are simply given by 
qW _ (_ T . \-in(») _ Qi 
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From Eqs. (2.33)-(2.35) it is evident that to give the frame components of the source we 
need to know only the perturbations 5p, 5p, and V % of the fluid. No metric perturbations 
are needed — in contrast to the coordinate components 5T® in which h 0i enters. This is 
important for our understanding of Mach's principle. 

III. FIELD EQUATIONS 

We have seen that the accelerations and rotations of LIFs with respect to the COF are 
determined in a general gauge by hoa and hoi components of the perturbations. We shall now 
write the perturbed Einstein's equations for the FRW backgrounds in a general gauge 2 . We 
then shall see later which gauges will enable us to determine instantaneously perturbations 
hoo an d hoi (separately from hki) in terms of matter perturbations. 

A straightforward way to express the perturbations of Einstein's equations is in terms of 
a physical cosmic time t and some convenient spatial coordinates x l of the FRW background. 
However, there are advantages in using conformal time 77, given by a(r)) drj = dt. Both t and 
i] are common in the literature and we shall thus give explicitly the perturbation equations 
in two forms — with t in Appendix A and with 77 in this section. 

A. Perturbed field equations with a conformal time r\ 

In terms of coordinates x^ = (r),x k ) the metric of the background is 

ds 2 = Jj^dx^dx" = a 2 e flI/ dx fl dx u = a 2 [drj 2 - f M dx k dx l ] , (3.1) 

where we introduced the conformally related static background metric e^ u by 

eoo = 1 , e ; = , e ki = -f k i(x l ) . (3.2) 

The components of a tensor, say W^, are related to those of the tensor W£ in x M = (t,x k ) 
coordinates as follows: 

Wg = W8, W? = a- l W?, W l = aW l Q , W{ = W l k . (3.3) 

2 Some of the equations presented here have been derived independently by Langlois (1994) in his PhD 
thesis and by P. Uzan and N. Dcruellc (private communication). The perturbed Ricci tensor components 
and the equations of motion have been written down by Bardeen (1980) [29] after a decomposition of 
the metric into scalar, vector, and tensor parts using his specific amplitudes. 



21 



Defining the dimensionless "relative Hubble constant" hy H = = — — a — aH, we can 

° J a dr\ a 

write the non-vanishing background Christoffel symbols as 

r _ u -pO U f "pm U X m r m "T m /o a\ 

where f ^ is given in Appendix A. The prime hereafter denotes the derivative with respect 
to rj. (Later, it will also be used to denote a coordinate change, but no confusion should 
arise.) The nonzero components of the background Einstein equations become 

C7o =Gl = - 2 {k + H 2 ) =Kp + K, G™ = G™ = \8™(k + H 2 + 2H') = -(Kp-A)5™ . 

CL CL 

(3.5) 

The linearly perturbed Einstein equations will be written in terms of the dimensionless 
perturbations h^ u of e^ u , 

ds 2 = {g^u + Sg^dx^dx" = a 2 (e^ + dx^dx v = 
= a 2 (1 + h m )drf + 2h 0k dr]dx k - (f kl - h kl )dx k dx l , (3.6) 

which means that Sg^ u = a 2 h^ u . It is important to emphasize that, in contrast to tensors 
like in Eq. (3.3), h^s are not components (in coordinates x M ) of used in (1.1) and 
Appendix A; as seen from Eq. (3.6), they represent the perturbations of the conformal 
static metric e^ u , whereas h^s represent perturbations of the physical background metric 
g^y. In a (l+3)-decomposition, i.e., in quantities h 00 , h i, h k i, we do not raise the index 
and we raise the spatial indices only with f kl ; thus h™ = f ml h i , h mn = f mk f nl h M , etc. 
The explicit relations between h^, W , or h^ v and h^, h" or are given in Appendix A. 

The perturbations of Einstein's equations in terms of h^ v can be obtained from equations 
in Appendix A. We shall introduce two special symbols which not only simplify the equations 
but are also helpful in suggesting particularly useful gauge conditions. We set 

/C = \nh m + \ (hi) - V n K ■ (3.7) 

The second quantity we shall employ is defined with the traceless part of h l k , 

h l Tk = h[- U l k hl . (3.8) 

We set 

% = Vih l Tk . (3.9) 
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Now we give Einstein's perturbation equations, separating SG l Tk , the traceless part of 5G l k , 
from the trace 5G™ which we combine with SG^ for a reason to be seen below. Thus, recalling 
that V 2 = f kl Vki, we have the following dimensionless equations 



a 2 K5T% 



a 2 5G° = l -V 2 h n n + khl - 2WC - l -V k T k , 
a ^G° k = ^V% + kh k0 + ^W kl h l 

+ l^kic - l - (%)' , 



a 2 { 5G° - SGI 



V 2 h 



oo 



:]<i [ -H j //,„, • - laK)' 
a I a 



(3.10) 



(3.11) 



(3.12) 



and 



<rn 1 5T l k - -slsr;; 



a 2 5G 



Tk 



V 2 h l Tk + kh l Tk + 



Tk 



2a 2 



a 2 h 



l Tkj 



+ f 



Irn 



V( m 7^) — -f mk S7 n T r 



1 

~a~ 2 
1 



^ jlm 
-.2 



a ( V( m /i fc )o - -fmk^nh 



+ -f lm (v mfe -i/ mfe V 2 J Uoo-^n 



(3.13) 



For completeness we also write down the equation which can be derived from Eq. (3.11): 



a 2 n 5T k 



1. 

2 



=:V 2 h k 



k-2a(-n)' 



hi 



l -V k V l h l - 2 -V k lC+ l -{T k )'. 



(3.14) 



This equation follows from Eq. (3.11) by using the relation 



5T° k = -f kl 6T* 



-k + a 



n 



hok- 



(3.15) 



In the case of perfect-fluid perturbations we define the local coordinate velocity by 



V k = 



dx k (r]) 
drj 



(3.16) 



Notice that V k is not equal to V k [defined in Eq. (2.36)] expressed in coordinates x M because 
V k is defined with respect to the conformal time. Since x k = x k , in both coordinates we 
have simple relations, 



V k = aV k , 

V = f V m = f aV m = -a~ x a V m = -a^V 

v n — Jnm v — Jnrn u 'V — u ynm v — u v n- 



(3.17) 
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Nevertheless, the fluid's 4- velocity components (2.37) transform as a general tensor like (3.3). 
In the case of perfect fluid the energy-momentum tensor perturbations become 

a 2 K 5T ° = a 2 K 5p , cl 2 k 5T° = 2 [k + H 2 - ft') (-V k + h k ^j , 

a 2 nST k = 2{k + H 2 -H')V k , a 2 K5T l k = -a 2 K5 l k 5p, (3.18) 

so that the left-hand sides of equations (3.12) and (3.13) are 

a 2 K (67$ - ST^j = a 2 K (Sp + 35p) , a 2 K (sT l k - U k 5T^ = . (3.19) 

We combined Einstein's equations in such a way that Eqs. (3.10) and (3.12) contain 
scalars under the transformation of spatial coordinates, whereas Eq. (3.13) involves tensorial 
quantities only. In the perfect-fluid case, the 'source' in Eq. (3.13) vanishes so that this 
equation represents propagation of a free gravitational field, i.e., of gravitational waves 
described by traceless quantities h l Tk . The first and the third term on the r.h.s. of Eq. (3.13) 
combine into a d'Alembert wave operator modified by the time dependence of the expansion 
factor a(rj). More generally, however, the perturbed fluid could be an imperfect fluid which 
includes shear viscosity. This can be described by an additional term in ST" given by a 
symmetric shear tensor <5E^ which is traceless, <5£^ = 0, and purely spatial in the fluid rest 
frame, WST,^ = (see, e.g., [26]). Then the shear would appear as a source in Eq. (3.13). 

B. Bianchi identities and conservation laws 

The perturbed contracted Bianchi identities, V^G^ = for p = and p = k, imply 

\(a?6G%y + -(6& -6G$- -V k 8& k 
cr a a 

^(a 4 5G° k y + V m <5<5™ - \n{p + p)V k h 00 = 0. (3.21) 
Replacing SG^ by 5T£ from the field equations we get the conservation laws for ST£: 

(5f^ + -(35f -5f^--V k 5f° k 

CL CL 

+ ^(P + P) (ahoo-^lCj =0, (3.22) 
4(a 4 5T°)" + V m 8T™ -\{p + P) V k h 00 = 0. (3.23) 
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In the case of a perfect fluid the conservation laws become 



{Sp)' + -(Sp + Sp) + -(p + p)V k (V k - h k Q ) 
a a 



0, 



1 



a 4 (p + p)(V k - h k0 ) + V k Sp + - (p + p)V k h 00 = 0. 



(3.24) 



(3.25) 



The first equation expresses the conservation of the mass-energy Sp. The second is the 
equation of motion; when the time-derivative term is negligible it represents the equilibrium 
condition between the gradients of pressure and gravitational potential |/ioo, which would 
be much harder to see in the synchronous gauge with /i o — 0. Until now, all equations 
have been in an arbitrary gauge. The next section is devoted to the choice of 'appropriate 
gauges'. 



IV. GAUGES 



A change of the gauge can be regarded as an infinitesimal coordinate transformation 
M + (^(x). Under the gauge transformations, the metric changes by the Lie 
derivative (e.g. [23]) as Ag^ = £ c g fXl/ = lim c ^ [g^(x') - g'^(x')] = C^ + C^u- The explicit 
formulas are given in Appendix B. 

Since gauge transformations contain four arbitrary functions, we can impose four gauge 
conditions. Regarding the field equations (3.10) - (3.13) we see instantly that four gauge 
conditions /C = = % decouple the first three equations from the rest. Equations (3.10), 
(3.11), and (3.12) determine directly the metric components hoo,hko, and /i™ from the in- 
stantaneous distribution of sources given by 5T$ , ST® and STq — ST™; no time integration 
is needed. Accelerations and rotations of local inertial frames follow then from (2.20) and 
(2.25). Such an instantaneous determination of local inertial frames is also possible by 
employing other gauges. We call these gauges Machian. 

The purpose of this section is to motivate and describe geometrically several Machian 
gauges, and to clarify what is the residual gauge freedom these gauges admit. For a compar- 
ison we shall also consider two typically non-Machian gauges — the synchronous gauge and 
the generalized Lorenz-de Donder, or "harmonic" gauge. In the next section these gauges 
will be used to analyze the field equations and the way they can be solved to determine 
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local inertial frames. In the Machian gauges we shall always restrict the spatial part of the 
metric by requiring the three gauge conditions T k = 0, where T k is given by Eq. (3.9). These 
conditions will be motivated first. 

A. Gauge conditions on the spatial metric 

We start beyond the linear perturbation theory. Smarr and York (1978) [48], in treating 
full general relativity as an evolution from given initial Cauchy data on a spacelike slice to 
the next slice, studied the kinematics of the observers threading the slices. The evolution is 
represented in terms of coordinates attached to these "coordinate observers" . Kinematical 
and dynamical effects can be suitably separated if a relative velocity of the coordinate 
observers, with respect to the (Eulerian) observers whose worldlines are perpendicular to 
the given slicing, is such that the shear of coordinates arising if one goes from one slice to 
the next is minimized. For a given slicing the relative velocity is determined by the shift 
vector and, therefore, Smarr and York require "the minimal-distortion" shift vector. This 
condition is equivalent to the equation 

D% = 0, (4.1) 

where %j = (det7)~ 3 7i J ' is the conformal 3- metric on a given slice, the dot denotes the time 
derivative, and denotes the covariant derivative with respect to the spatial metric 7^ 
induced on given slicing by four-dimensional metric g^. The condition (4.1) is a natural 
choice from a number of points of view. We refer the reader to the original paper [48] for 
the details; here we wish to make just a few comments. In the weak-field limit in the wave 
zone, condition (4.1) generalizes and includes the well-known "transverse-traceless" gauges 
of Arnowitt, Deser and Misner (1962) and of Dirac (1959) (see, e.g., [1]). In the linearized 
gravity in generally curved coordinates with = g^+h^, where goo = 1, goi = 0, d t g~ij = 0, 
the condition (4.1) implies d t D^h,T%j = 0, where tiTij = — \ hg~ij, h = g^hij, and is 
the spatial covariant derivative with respect to g~ij. This is analogous to the radiation (or 
Coulomb) gauge condition in electrodynamics. In stationary spacetimes with a timelike 
Killing vector £ a the gauge condition (4.1) is satisfied if the slicing is carried into itself by 
the £ a isometry and £ a is tangent to coordinate observers. It is interesting to consider more 
general slicings. Choosing in Schwarzschild spacetime the slices orthogonal to the geodesies 
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of particles freely falling from rest at infinity, one finds that the condition (4.1) implies 

ds 2 = (1 - 2M/r)dr 2 - 2(2M/r)*drdr - dr 2 - r 2 [dQ 2 + sin 2 9ip 2 ), (4.2) 

i.e., one obtains the time-independent form (4.2) with the spatially flat metric on the slices. 3 
Now in our case of the perturbed FRW metric it is easy to see that the conformal 3-metric 
is 7^ = f~*[fij — h Ti j], f = det(fij), and hnj is given by (3.8). Therefore, the condition 
(4.1) implies T k = 0, with T k given by (3.9). It is the last condition which converts Eq. 
(3.11) into the equation for h k o without the terms depending on the traceless part of h k \. 

Motivated by the analysis above, we shall assume that, in fact, a slightly stronger condi- 
tion, 

% = V$ Tk = 0, (4.3) 

is satisfied. This just means that the spatial coordinates are restricted on an initial slice and 
this restriction is then maintained by the original condition. Notice that (4.3) is covariant 
under 3-dimensional coordinate transformations within chosen slices. 

B. Gauge conditions on the time slicing 

The three gauge conditions T k — 0, k — 1, 2, 3, do not restrict the time coordinate, i.e., 
the slicing by spatial hypersurfaces x° = constant. We thus supplement them with the 
fourth gauge condition fixing the slices. In order to understand its geometrical meaning, 
we first calculate the geometrical quantities characterizing the slices. Using the perturbed 
FRW metrics in a general gauge, we find that the unit timelike vector field n M orthogonal 
to each slice is given by 

W = + 5W, (4.4) 

where n M = (a -1 , 0, 0, 0) , 5n^ = a -1 ^— |/ioo ; ^o) • Calculating the expansion 9 = of the 
congruence of timelike curves that meet the slices orthogonally, we find 

3 More recently, (4.2) was rediscovered as a technically suitable form of the Schwarzschild metric for de- 
scribing the Hamiltonian dynamics for spherically symmetric gravitating shells by Kraus and Wilczek 
(1995) [49] without any geometrical argumentation. 
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where /C is given by Eq. (3.7), and 

V = (4.5) 
For the shear of the congruence we obtain a a p = Sa a/ 3, 5a 00 = 0, 5a 0i = 0, 

Sa i:j = -aV(Jij)o + ^afijVrJi™ + ^ah' Tij 

~ 1 1 ~ 

= -aV(ih j)0 + -afijV + -ah' Tij . (4.6) 

The uniform- Hubble- expansion gauge, introduced by Bardeen (1980) [29], but apparently 
not used much later (though see [50]), requires 86 = 0, i.e. 

)C = ~h 00 + 1 -(h:y-V = 0. (4.7) 

This gauge condition is again motivated by the popular choice of the 'constant mean curva- 
ture slices' in the full theory (the trace of the extrinsic curvature tensor - the mean curvature 
- of a spacelike hypersurface with normal is K = —n? ). The condition (4.7) thus means 
that we choose such a time coordinate in the perturbed FRW universes that the extrinsic 
curvature of the rj = constant hypersurfaces is the same as in the unperturbed universe, 
i.e., it is constant along each hypersurface. Much work has been done on the existence and 
properties of such foliations (see [51] for the recent review and references). 

The gauge condition (4.7) for time slicing combined with the gauge conditions (4.3) for 
spatial part of the metric will be called the Mach 1 gauge. We have not found it in the 
literature, although the gauge conditions (4.3) and (4.7) were used separately. 

Another basic geometrical object associated with a spacelike slice is its intrinsic curvature, 
the simplest measure of which is the intrinsic (3-dimensional) scalar curvature. In the 
perturbed FRW universes K = K + 5K, H = -g, and 

STZ = "^(V 2 ^ + 3kh n n ) + ^V„T" (4.8) 

When the gauge condition (4.3) is combined with "the uniform-intrinsic-scalar curvature" 
condition 51Z = 0, i.e. 

V 2 K + Zkhl = 0, (4.9) 
we speak about the Mach 2 gauge. A stronger version - the Mach 2* gauge - requires 

K = 0. (4.10) 
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Another possible condition for the choice of slicing is ViVjK lj = 0, where K tj is the extrinsic 
curvature tensor. Nothing appears to be known about this choice in the nonlinear context. 
In our formalism this condition reads [using Eq. (4.6) for the shear] 

= V l V j 5a l3 = -\a{V 2 + 3k)Vih l + ^aV l T/, (4.11) 

which justifies the name "minimal-shear hypersurface condition" suggested by Bardeen [29]. 
Combined with the conditions (4.3) the last equation implies the gauge condition 

(V 2 + 3k)Vih l = 0. (4.12) 

The Mach 3 gauge is defined by the gauge conditions (4.3) and (4.12). Its stronger version, 
the Mach 3 * gauge, 

V = V l h l = 0, (4.13) 

combined with (4.3), has been called the Poisson gauge by Bertschinger [26] in 1995. He 
analyzed its advantages for physical interpretation of cosmological perturbations, in partic- 
ular, as compared with the synchronous gauge. The same gauge has already been proposed 
in 1994 by Bombelli, Couch, and Torrence [52] who called it the "cosmological gauge" . 
We also mention the standard synchronous gauge, still used most commonly in cosmology, 

h 00 = h 0i = 0, (4.14) 

and, in more detail, the generalized Lorenz-de Bonder gauge (frequently also called the 
harmonic gauge - cf., e.g., [3], recently [53]), 

%S(y/=gg^) = 0. (4.15) 

This has been extensively used in a number of problems, in particular, in weak-field approx- 
imations dealing with equations of motion and gravitational radiation (see, e.g., [54]), but 
not in cosmology. The gauge conditions (4.15) for /i = imply 

-^i^ + l{h 00 + h n n )' + H{3h 00 + h n n ) =0, (4.16) 

for /i — k we get 

"Vip + \^k(k - 3M + h' ok + Mhok = 0. (4.17) 
We shall now analyze the residual gauge freedom which the gauges introduced above admit. 
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C. Gauge-fixing and residual gauge freedom 



From relations (B3) - (B6) in Appendix B we readily obtain the changes of the geometrical 
quantities defined in (3.7), (3.9), (4.5), and (4.6) under gauge transformations: 

AT k = A(V t h l Tk ) = -(V 2 G + 2k( k + ^V fc V,C'), (4.18) 
AV = A(V^) = -V 2 C° - aV k C\ (4.19) 

(X 



A/C = -- 

a 



V 2 C° + 3a 2 (^l C 



(4.20) 



ATI = -1(V 2 + 3A;)(V^ + 3-C°), (4.21) 
6a z a 

A{V k V l 5a kl ) = -jj(V 2 + 3A;)V 2 C°. (4.22) 

We shall discuss first the minimal- distortion spatial gauge condition: T k = 0. 

Starting from a general gauge, we reach the required condition by purely spatial gauge 
transformations given by ( k , which satisfy the inhomogeneous equation with given l.h.s. 
AT,. The residual gauge freedom is determined by ( k solving the homogeneous equation 

V 2 C fe + 2k( k + 1 V fe V,C' = 0. (4.23) 

There are solutions of this equation given by linear combinations (with time- dependent 
coefficients) of the conformal Killing vectors in the constant-curvature spaces S 3 ,R 3 ,H 3 . To 
see this, recall that in 3-dimensional space a conformal Killing vector satisfies 

V,C* + V fc Ci = pkiV n C. (4.24) 

Since spaces of constant curvature are conformally flat, they admit ten linearly independent 
conformal Killing vectors as E 3 (see, e.g., [55]). Their explicit forms are given in Appendix C, 
where their relationship to the scalar and vector hyperspherical harmonics is also elucidated. 
Among the ten conformal Killing vectors, six are pure Killing vectors, ^\ A \A = 1,2, ...,6; 
the remaining four ip\ ,B = 1, ...,4 do not reduce to the Killing vectors. In spaces of 
nonvanishing constant curvature, can be written as gradients of scalars: 

fc = ±l: 4 B) =d t Q B = V t Q B , B = l, ...4. (4.25) 
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The four scalar fields, Q B , are equal, up to a multiplicative constant, to the following four 
scalar harmonics (see Appendix C for details): 

k = +1 : Q( L =i, 1=0, m=o) - cosx, Q(L=i, i=i, m=-i, o,+i) - sin xY"i m (6>, V 7 ) , ( 4 -26) 
fc = -1 : Q(A=2t, i=0, m=0) - COshx, Q(A=2t,J=l,m=-l ) 0,+l) - sinh x *im(0» ¥>)■ (4-27) 

Taking the divergence of Eq. (4.24) and using (A5) and (A6) in Appendix A to commute the 
derivatives, we obtain Eq. (4.23). Hence, any conformal Killing vector ( k solves Eq. (4.23). 
In open universes, all such solutions (k diverge at infinity (x — > oo), except for translations 
in flat (k = 0) universe when (k are constant in Cartesian-like coordinates. We now prove 
that there exist no bounded solutions of Eq. (4.23) other than conformal Killing vectors in 
S 3 and translations in E 3 . To prove this we decompose Cfc into a gradient of a scalar and a 
transverse vector: 

Cfc = V fc Z + CTfc, V fc C™ = 0. (4.28) 

This decomposition is unique up to Z — > Z+constant, if for open universes we require (k t° 
decay asymptotically so that / V k (,kdV converges [23]. Substituting then (4.28) into (4.24) 
and commuting the derivatives, we find 

V 2 V 2 Z + 3W 2 Z = 0. (4.29) 

In S* 3 the only smooth solutions of the equation 

V 2 Q + 3Q = (4.30) 

are given by the linear combination of the four scalar harmonics (4.26), the gradients of 
which give the conformal Killing vectors. In closed space the solution of (4.29) is thus, with 
(3 sit) arbitrary, 

4 

Z = J2Mt)Q B + Z (t). (4.31) 

B=l 

In H 3 equation (4.29) becomes 

V 2 Q - 3Q = 0. (4.32) 

The solutions are four scalar harmonics given in (4.27) - these, however, diverge at infinity. 
In open universes the only solution of (4.32) leading to asymptotically well-behaved Cfc is 
Z = Zo(t), the gradient of which does not contribute to Cfe- Substituting now for in 
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Eq. (4.23) the decomposition (4.28), and regarding the above results for Z, we find that 
Eq. (4.23) reduces to the equation V 2 C,Tk + 2/cCt/c = 0. In open universes this equation 
does not admit any asymptotically well-behaved solutions, except for ( Tk = constant for 
k — 0. In a closed universe the equation is equivalent to the Killing equation (Appendix 
C). Hence, under the assumption that the vector is bounded, the condition (4.3) fixes 
the spatial coordinates uniquely in H 3 ; and in E 3 the remaining gauge freedom is just 
= Ea=i fA(t)Ctf''\ corresponding to a time- dependent linear combination of translations. 
In S* 3 the residual gauge freedom is given by a linear combination of ten conformal Killing 
vectors (six Killing and four conformal Killing): 

c = J2 Mm {A)i + £ M*w (B)i - (4.33) 

A=l B=l 

We now discuss the three Machian gauges successively. 

1. Mach 1: Uniform- Hubble- expansion gauge 
From Eq. (4.20) we see that the residual freedom in (° is given by the solutions of 

V 2 C° + 3a 2 (^j C° = 0. (4.34) 

Multiplying Eq. (4.34) by (° and integrating by parts over a domain V, we find 

/ c°v 2 c°rf (3 V = / (°v k (°ds k - / f^^v^v 

Jv JdV Jv 

= - J 3a 2 (^j (C°)V 3 V, (4.35) 

where S 3) V = d 3 x, f = det(f kl ); a 3 S 3) V, a = a(x°) is the proper volume in a slice x = 
constant. Taking V to be all space, then the integral over the boundary vanishes in open 
spaces because of the boundary condition on (°, and it is zero in closed spaces because there 
is no boundary. Therefore, 



I 



f'VkCVtCd^V = 3a 2 J ((°) 2 d^V. (4.36) 



IV 

The factor on the r.h.s. can be rewritten using the FRW background equations: 
A{t)=3a 2 (^\ =3a 2 H = 3k-\a 2 K(p + f) = -3a 2 {H 2 -\X)-\a 2 K(p + 3f). (4.37) 



In all standard models the strong energy condition p + 3p > is valid so that A < 
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(H 2 — | A > is satisfied in realistic models). In inflationary universe models with p + p 



(A = 0), the function A(t) < for open universes. In all these cases the r.h.s. of Eq. 
(4.37) is nonpositive, whereas the l.h.s. is non-negative. Therefore, the only solution of 
Eq. (4.34) is (° = 0. In the standard inflationary model with k = 0, p + p = 0, we have 
.4 = 0, and (° = (°{t) is an admissible solution of Eq. (4.34) which is bounded and has a 
vanishing gradient (reflecting the higher symmetries of de Sitter space to which the FRW 
models reduce). If k — +1 and p + p = 0, the relation (4.37) turns Eq. (4.34) (for any A) 
into V 2 C° + 3C° = 0, which is Eq. (4.30), the solutions thus being 

4 

(° = J2Mt)Q B , (4-38) 

B=l 

where 4 scalar harmonics Q B are given in (4.26), Ob are arbitrary. 

Let us summarize. Assuming bounded, the Mach 1 gauge fixes the coordinates 
uniquely in the open universes with k — — 1, and for k = it determines the spatial 
coordinates up to time-dependent translations, ( k (t), whereas the time slicing is unique if 
the background matter satisfies the strong energy condition; x° can be shifted by (°(t) m 
the inflationary universe. In closed universes the spatial coordinates are determined up to 
the time-dependent motions (4.33) given by the Killing and conformal Killing vectors; the 
time slicing is unique in the standard backgrounds with the strong energy condition satis- 
fied. In the inflationary backgrounds the time can be shifted by (° determined by Eq. (4.38). 

2. Mach 2: Uniform- scalar- curvature gauge 
Requiring the scalar 3-curvature of the time slices to be equal to the background values fixes 
the gauge up to the transformations satisfying [see Eq. (4.21)] 

(V 2 + 3A0(V ; C' + 3-C°) = 0. (4.39) 
a 

Assuming again the condition (4.3), we restricted ( l already by Eq. (4.23), which implies 
the divergence V;C' to satisfy (V 2 + 3k)Vi( l = 0. Equation (4.39) thus reduces to (assuming 
a 7^ 0) (V 2 + 3k)(° = 0. As discussed above [cf. (4.30) or (4.32)], the only bounded solutions 
are (° = if k = —1, (° = (°{t) if k = 0, and (° is given in terms of Q B for k = +1. 

Mach 2*: The traceless gauge. - The gauge condition ft,™ = implies the previous one, 
and is stronger. Indeed, regarding Eq. (B6), we see that the residual gauge freedom is given 
by VfcC fc + 3^C° — 0> which determines (° in terms of ( k (assuming a^O). With the gauge 
conditions (4.3), (° = in open spaces, in S 3 the residual freedom in ( k is given by Eq. 
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(4.33) which implies V k ( k = Eb=iM*) v2( 2 B [see Eq. (4.25)], and thus leads to [using 
(4.30)] 

4 

(° = (a/a)Y,Mt)Q B , (4.40) 

B=l 

(3 Bif) are arbitrary. 

3. Mach 3: The minimal-shear gauge 
As seen from Eq. (4.22) this gauge condition allows the transformations restricted by 

(V 2 + 3ife) ( V 2 C° - a 2 V k ( k ) = 0, (4.41) 

which, using; C fc that satisfy (4.3), reduces to (V 2 + 3k) V 2 C° = 0. This is the same as (4.29). 
In open spaces the only bounded solutions are (° = C°(t). In closed spaces 

4 

C° = $> B (t)Q B + a (t), (4.42) 

B=l 

where cr's are arbitrary, and Q B given by Eq. (4.26). 

Mach 3*: The Poisson gauge - The condition Vih l = admits a smaller freedom. Equa- 
tion (4.41) becomes 

V 2 C° - a 2 V k ( k = 0, (4.43) 
which, after substituting for ( k from Eq. (4.33), for k — +1 gives 

4 

(° = Y,a 2 Mt)Q B + a°it), (4.44) 

B=l 

where a is arbitrary but (3 sit) are the same functions as those in ( k in Eq. (4.33) - in 
contrast to Eq. (4.42) where cr's are independent of /3's. In open universes the only residual 
freedom in the choice of time is given by arbitrary C°(0- 

Regarding the gauge freedom in T k = 0, we see that the Poisson gauge in the case k = — 1 
fixes the spatial coordinates uniquely; the time coordinate is fixed up to (°(t). In the case 
k — the spatial coordinates are fixed up to translations C(t) and time shifts C°(0- ^ n ^ ne 
closed case the freedom in spatial coordinates is determined by linear combinations of the 
Killing and conformal Killing vectors (4.33), whereas the time coordinate by the combination 
(4.44) of scalar harmonics Q B . Hence, in the closed case there are 11 arbitrary functions of 
time which represent the gauge freedom. 

34 



These results are at variance with Bertschinger's statement [26] that there is "an almost 
unique transformation from an arbitrary gauge to the Poisson gauge" . Clearly, Bertschinger 
does not consider the possibility that his (3 solves equation (V 2 + 3k)V 2 (3 = [i.e. our 
Eq. (4.29)] which preserves the gauge condition Vih l Tk = 0. Solutions for f3 for k = ±1 are 
as in Eq. (4.31), where Z (t) indeed has no effect but the terms containing Q B do have an 
effect - not only on (k but also on (° as described above. 

Let us now mention the gauge freedom in two typical 'non-Machian' gauges. 

4- Synchronous gauge 

From Eqs. (B3) and (B4) it is immediately seen that the synchronous gauge admits the 
well-known residual freedom given by transformations satisfying 



C° = 0, V,C° = a%, 



(4.45) 



which imply 



C° = CV), 



dt 



v,cV) + ^ 1 )- 



(4.46) 



Functions (°( xt ) an d %i(x l ) are arbitrary. The gauge freedom is the same for both open 
and closed universes. 



5. The generalized Lorenz-de Donder gauge 
Requiring the gauge conditions (4.16) and (4.17) to be satisfied, we can use relations (B3) - 
(B6) to find out the residual freedom in this gauge. It turns out to be restricted by 

V 2 C° - a 2 C° - 3aaC° + 3(aa + a 2 )(° - 2aaV l ( l = 0, (4.47) 

V 2 ^ - a% + 2kC k - haaCu + 2(d/a) V fe C° = 0. (4.48) 

The only feasible way to solve this coupled system appears to be the use of harmonics, but 
here we shall just restrict ourselves to noticing that, for a slowly changing expansion factor 
(d, a small), the system turns just into two decoupled wave equations, 

lv 2 C° - C° = 0, - 2 V\ k - C, + = 0. (4.49) 
a 2 a 1 a 1 

In the flat case these are just wave equations in flat space with coordinates ax 1 (which 
give the proper lengths in k = universes). The gauge freedom is thus analogous to the 
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freedom of the Lorenz gauge in electrodynamics. Any solution of a wave equation can be 
characterized by its Cauchy values - here C ^*); (°( xt ), (k(x l ), and Cfc(x*), i.e., by eight 
functions of spatial coordinates. 

Summarizing, we find that the Machian gauges are substantially more restrictive than 
the synchronous gauge and the generalized Lorenz-de Donder gauge. The last two gauges 
admit transformations characterized by several (two and eight) arbitrary functions of three 
variables - of the spatial coordinates x\ All the Machian gauges admit only several arbitrary 
functions of time. In some cases they fix the coordinates uniquely. An arbitrary additive 
function of time, C°(0> represents just the changes of the units of time: dt' = (1 + (°)dt. 
The spatially homogeneous changes of x 1 by Q{t) describe just the shifts of the origin of 
spatial coordinates. Otherwise, all three Machian gauges fix the coordinates uniquely in the 
hyperbolic universes H 3 as a consequence of boundary conditions at infinity. The remaining 
functions of time in spherical universes S* 3 will be interpreted in the following. 

D. Integral gauge conditions in closed universes with standard spherical topology 

In the closed spherical spaces the spatial coordinates x % are fixed up to the transformations 
x % — > x n = x l + C where ( l is given by a linear combination of six Killing and four (proper) 
conformal Killing vectors of S 3 , in which the coefficients are arbitrary functions of time. 

In order to acquire an insight into the effects such coordinate changes can produce, con- 
sider an unperturbed FRW universe with standard spherical topology (k = +1). Transform 
the metric in the hyperspherical coordinates (1.3) by a gauge transformation generated by 
one translational, one rotational, and one conformal Killing vector which have the simplest 
forms in the hyperspherical coordinates: Q T = (cos#, — cot x sin 0, 0), Q ot = (0, 0, l),Cconf = 
(— sin x, 0,0) (cf. Appendix C). Admitting the time-dependent coefficients, the transforma- 
tion has the form 

x' = X + a(t) cos0 - 7(t) sinx, 

0' = 0-a(*)cotxsin0, <p' = (p + P(t), (4.50) 

which can easily be inverted since a, (3, 7 are small. In addition to the transformation (4.50) 
we consider a change of the time coordinate (time slicing) allowed by our Machian gauge 
conditions in closed universes. 
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Hence, we take C° of the form (4.42) because other possibilities (4.38), (4.40), and (4.44) 
are included in (4.42). However, since in (4.50) only the simplest conformal Killing vector 
enters, it is sufficient to take only those time transformations which are associated with this 
vector and with the shift of the time origin which is also allowed by Eq. (4.41): 

t' = t + a 2 [6(t) cos x + cr(t)] , (4.51) 

where for convenience the expansion factor is pulled out. Under the transformations (4.50) 
and (4.51) the standard FRW metric with k — +1 becomes 

ds 2 = [1 - (2a 2 5Y cos X ' - (2aV)'] dt' 2 

- a 2 [1 + 2(7 + add) cosx' + 2aaa] [dx' 2 + sin 2 x\d9' 2 + sin 2 9'dip' 2 )] 
+ 2a 2 {[dcos6>' + (5 - 7) sin*'] d*' 

- a sin x' cos x' sin 9'd6' + (3 sin 2 X ' sin 2 9'dip'}dt'. (4.52) 

Since a, a are, in general, time-dependent and h' 0Q , h' 0i nonvanishing, the frames associated 
with x' 1 6' ■> V 9 ' fixed are noninertial in general, and the inertial frames, momentarily at rest 
with respect to them, are seen to have the acceleration [cf. Eq. (2.20)] 

a* = -(l/a 2 ) { [2a 2 (5 - 7)] ' sin x ' + (a 2 «)' cos 9'} , 

a' e = -(l/a 2 )(a 2 «)-cotx'sin#', a ,lp = -(l/a 2 )(a 2 /3)', (4.53) 

and to rotate with the vorticity [cf. Eq. (2.25)] 

uj' x q = a 2 a sin 2 X ' sin 9', ui' xip = a 2 f3 sin X ' cos X ' sin 2 9', 

= - a 2 (3 sin 2 X ' sin 9' cos 9', (4.54) 

i.e., with the vorticity 3-vector (2.31) given by 

u' x = (l/a)/3cos9', uj' 6 = (l/a)$ cot x sin 9', J* = (l/a)a. (4.55) 

The above results are easily understood: time-dependent rotations in the tp direction, with 
X = 0, and the f axis fixed, imply nonvanishing accelerations in this direction only, whereas 
the corresponding vorticity vector has no (p component. The translational Killing vector 
which for, say, 9 = represents rotations in the X direction (with X = n/2 fixed), leads to 
accelerations only in the X direction, and the vorticity vector in the if direction. 
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None of these acceleration or vorticity vectors can be compensated by an allowed change 
(4.51) of time slicing. As expected, the shift of the time origin, a(t), does not enter Eqs. 
(4.53)-(4.55). However, the effect of the transformation generated by the conformal Killing 
vector, which appears only in the x component of the acceleration, can be annulled by 
choosing 8 — 7. Nevertheless, this condition does not remove the effect of both 5 and 7 in 
the conformal factor of the spatial background metric. Metric (4.52) implies a nonvanishing 
trace of the form (omitting "the time shift" a) 

K = ~K = 6 (7 + aaS) cos X '. (4.56) 

Therefore, the spatial metric differs from the canonical metric of a homogeneous and isotropic 
3-sphere. This metric is preserved only by transformations representing real symmetries, i.e., 
those generated by the Killing vectors. 4 

We have described the particular effects of the gauge freedom corresponding to the Killing 
and (proper) conformal Killing vectors in order to show their different character. In a gen- 
eral, linearly perturbed, FRW universe, the metric will be much more complicated than that 
of Eq. (4.52). If it contains terms appearing in (4.52) (and those corresponding to other 
Killing and conformal Killing vectors), they can, of course, be removed by gauge transfor- 
mations of the form (4.50) and (4.51). Since a natural goal in a relativistic perturbation 
theory is to fix the gauge at the end as uniquely as possible, we shall now require, in all 
Machian gauges, additional gauge conditions which exclude the freedom corresponding to 
the gauge transformations generated by (proper) conformal Killing vectors. However, we 
leave the freedom corresponding to the Killing vectors since these exhibit the symmetry of 
the background universe at any given time. 

The trace (4.56) is, at given time, proportional to cosx [we omit primes in metric (4.52)], 
i.e., just to the first of the scalar harmonics in Eq. (4.26). Neglecting the time shift ~ <r(t), 
the perturbation /i o = ^00 in (4-52) is also proportional to this harmonic. The term in 
h 0i (— ahoi) in (4.52) corresponding to the same conformal transformation of the spatial 
coordinates and of the time slicing is proportional to the gradient of cosx but the scalar 
V = Vih l is again proportional to this harmonic. 5 The metric perturbations which arise 

4 As in flat space, the canonical flat-space metric in Cartesian coordinates is preserved only by rigid trans- 
lations and rotations, but not by dilatations. 

5 We do not consider the scalar V/Vfc/i^since the gauge condition (4.3) guarantees that it vanishes in all 
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or may be removed by gauge transformations generated by the conformal Killing vectors 
(i.e., the vectors of the form ( l = Y1b=i Pb^^ 1 = Y1b=i Pb^iQ b ) will be eliminated by 
the following integral gauge conditions. These will be imposed at all times: 

/ h n n Q B d^V = [ h 00 Q B d^V = [ VQ B d^V = 0, (4.57) 
Js 3 Js 3 Js 3 

where h™,h 00 ,V = Vih l are functions of all spacetime coordinates, harmonics Q B (x,0,(p) 
are given in Eq. (4.26). The integral gauge conditions (4.57) require that spatial scalars 
h™, h 00 , and V are orthogonal to the 4-dimensional functional space spanned by Q B , i.e., by 
harmonics which are eigenfunctions with zero eigenvalues of the operator (V 2 + 3) in S 3 . 
In Section V we shall notice that conditions (4.57) are closely related to Traschen's integral 
constraints [37], [38] which restrict perturbations of energy-momentum tensors representing 
sources. In this way we make sure that conditions (4.57) do not restrict physics. 

E. Machian gauges in closed spherical universes: Summary 

After adopting the integral gauge conditions (4.57), the gauge freedom in all three 
Machian gauges becomes transparent and simple: It reflects the proper (Killing) symmetry 
of the background universe at any fixed time. The minimal-distortion shift gauge condition 
(4.3), together with integral gauge conditions (4.57), which we assume in all Machian gauges, 
fix the spatial coordinates uniquely up to transformations x % — > x h = x l + C with 

= X>a(^ (A) \ (4-58) 

A=l 

where are six spatial Killing vectors three (quasi)rotations and three 

(quasi)translations. 

In the Mach 1 gauge the time slicing is unique even without requiring integral gauge 
conditions (4.57) if the background matter satisfies the strong energy condition. After 
requiring (4.57), it is unique also in the inflationary backgrounds. In the uniform-scalar- 
curvature gauge (Mach 2), with (4.57) satisfied, the time coordinate is unique, the same 
being true for the special case of the traceless gauge (Mach 2*). Finally, in the minimal- 
shear gauge (Mach 3) and its special case of the Poisson gauge (Mach 3*), the adoption of the 

Machian gauges. 
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integral conditions (4.57) leaves the only freedom in Eqs. (4.42) and (4.44) to be (° = <J°(t), 
i.e., the time coordinate is fixed up to trivial, 'universal' shifts depending just on an arbitrary 
function of time. This, as noticed below Eqs. (4.53)- (4.55), does not influence accelerations 
and rotations of local inertial frames. Therefore, all our Machian gauges fix coordinates 
uniquely up to the 'time-dependent' Killing motions (4.58) of spatial coordinates. The 
Machian gauges with the integral gauge conditions (4.57) are thus determining coordinates 
both more restrictively and more plausibly than the synchronous and generalized Lorenz-de 
Donder gauge. 

F. On closed hyperbolic and flat universes 

Although in this work we generally assume the cosmological backgrounds with standard 
topologies only, and thus with geometries which are homogeneous and isotropic also globally, 
in this intermezzo we consider 3-dimensional backgrounds represented by closed flat (k = 0) 
and hyperbolic (k = —1) 3-manifolds. Finite universes with multiconnected topologies have 
become popular in recent years in the light of new theories extending general relativity, and 
in the view of a possibility (in principle) to determine the topology of our universe by means 
of cosmic microwave background observations, or from the distribution of distant sources. 
A comprehensive, nice review containing many references appeared recently [56]. 

Globally, these universes are different; in particular, they admit smaller families of con- 
tinuous symmetries. Closed hyperbolic 3-manifolds do not have smooth Killing vectors at 
all [56] and do not possess nontrivial solutions of the equation V 2 + 3k<p = 0. Indeed, 
multiplying this equation by 0, integrating by parts over a 3-dimensional domain V [cf. Eq. 
(4.35)], we get 

/ (pV^dSj - [ Vj0V J > d^V + 3k I <p 2 d^V = 0. (4.59) 
Jav Jv Jv 

Taking T> to be whole closed space, the first integral vanishes, because there is no boundary, 
and since both the second and the third integrals are non-negative, Eq. (4.59) for k = — 1 
can be satisfied only with = 0. Hence, Eq. (4.32) has only solutions Q = 0, so Z in Eq. 
(4.28) does not contribute to ( k . Analogously, equation for the transverse part of ( k , after 
being multiplied by (Tml km , and integrated by parts, becomes 

/ l mk (Tm^( Tk dS 3 - [ 1 mk Vj( Tm ^(Tkd {3) V + 2k [ ^(rUrkd^V = 0. (4.60) 

JdV JV JV 
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Again, taking V to be all closed space, the first term is zero, and as the other integrals are 
non-negative, (4.60) for k — — 1 is solved only by ( Tk = 0. It is easily seen that, except 
for trivial shifts a(t) in the Mach 3 gauge, in closed hyperbolic universes our instantaneous 
Machian local gauge conditions fix coordinates uniquely, without integral gauge conditions 
being imposed. 

For closed flat (k = 0) universes, Eq. (4.59) implies = 4>(i) and Eq. (4.60) gives 
VjCrm = so that C = Xm=i /^(OCt^) where ^ t ^\A = 1,2,3 are translation Killing 
vectors. Hence the condition (4.3) determines the spatial coordinates up to arbitrary time- 
dependent linear combinations of translations. That the rotational Killing vectors are glob- 
ally ruled out can be well understood in the simplest example of a compact flat 3-manifold 
- a 3-torus T 3 . All closed flat 3-manifolds are given in Fig. 26 in [56]. In fact, only the 
3-torus admits globally a 3-parameter family of translational symmetries given in Cartesian 
coordinates by the 3 independent constant Killing vectors. Considering just T 3 , we find that 
the gauge condition (4.3) fixes spatial coordinates up to 

3 

x i ^ x n = x i + C(t ) + £ f A (t)^f } \ (4.61) 

A=l 

where £ t ^\ A = 1, 2, 3, are 3 translation Killing vectors. In the Cartesian-type coordinates 
these can be chosen as = (1, 0, 0), etc., so that the transformation generated by them is 

x = x'-f {1) (t), y = y'-f ( 2 ) (t), z = z'-f {3) (t). (4.62) 

This brings the FRW background metric with k = into the form 

ds 2 = dt 2 - a 2 {dx' 2 + dy' 2 + dz' 2 ) + 2a 2 {f (1) dtdx' + f {2 )dtdy' + f {3) dtdz'). (4.63) 

The acceleration of the local inertial frames with respect to the frame with x', y', z' fixed is 
thus given by [cf. Eq. (2.20)] 

a' x = -(l/a 2 )(a 2 /(i))', ^ = - (1 / a 2 ) (a 2 f m ) , a' z = -(l/a 2 )(a 2 / (3) )'. (4.64) 
Since the gradients h' ok l are vanishing, local inertial frames do not rotate. 



41 



V. FIELD EQUATIONS, INTEGRAL CONSTRAINTS, SOLUTIONS 
AND INERTIAL FRAMES. 



We now turn to the equations for perturbations in the Machian gauges. Here we pay 
attention to the Mach 1 gauge. Its choice of the constant mean curvature slices is most 
natural from the perspective of the full nonlinear theory. Moreover, the structure of the field 
equations for linear perturbations and their solutions do not differ significantly for the Mach 
gauges considered. The equations in Mach 2 and 3 gauges and in the generalized Lorenz-de 
Donder gauge are briefly discussed in Appendix D. Whenever solutions are known in terms 
of Green's functions, we write them down. They can be used to determine the accelerations 
and rotations of local inertial frames. Alternatively, solutions in terms of harmonics [22], 
[57] can be obtained by direct calculations, but they will not be studied in the present work. 



A. Field equations in Mach 1 gauge 



The minimal-distortion shift condition (4.3) is combined with the constant mean external 
curvature condition, i.e., with [c.f. Eq. (4.7)] 



^d/ioo + \oh n n -V = 0, V = Vih l . 



(5.1) 



In addition, we impose integral gauge conditions (4.57). As a consequence of the differential 
gauge conditions (4.3) and (5.1), the field equations (3.10), (3.11), (3.12), (3.13), and (3.14) 
simplify considerably: 

W 2 h n n + Zkh n n = 3a 2 /€(5f °, 



V 2 h k0 + 2kh kQ + -V k V = 2cl 2 k5T°, 



W 2 h m + 3a 2 ( - ) h 00 = cl 2 k(5T% - 8T%). 



Instead of Eq. (5.3) we may, equivalently, consider the equation 



V 2 h k - 2 



k-2a[-n 

a 



h k + -V k V = -2a 2 K5f^ 
3 



(5.2) 
(5.3) 

(5.4) 
(5.5) 



in which STq plays the role of a source. The field equation (3.13) can be written in the form 
V 2 h l Tk - 2kh l Tk - i (a 3 /^)' + F{h 00 ,K,h ok: h ok } = -2a 2 K5T l Tk , (5.6) 
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where J- {■■■'} denotes terms linear in the quantities in the brackets and in their spatial 
derivatives. If T and STip k are known, the last equation is a wave-type equation for h l Tk . 

It is remarkable that in this gauge neither of equations (5.2)-(5.5) contain any time 
derivative. All four of these equations are elliptic equations for h™, h k Q, and /iqo when the 
right-hand sides are given. The first two are standard constraint equations, and the third 
became elliptic equation for h 00 as a consequence of the gauge conditions. Another remark- 
able feature of (5.2)-(5.5) is that, with ST®, ST®, ST™, STq given, they represent a completely 
separated system of four equations for, subsequently, h™, h k o, and hoc,. V = Vih l in Eq. (5.3) 
is governed by a separate equation. Applying V fc on Eq. (3.11) and commuting derivatives, 
we obtain 

V 2 V + ZkV + V 2 /C = -a 2 KV k ST®, (5.7) 

2 

which, in gauges for which /C = 0, turns into 

V 2 V + ZkV = ^-a 2 nV k ST®. (5.8) 

This has exactly the same form as (5.2) for h™. With ST® given, we can solve (5.8) for V 
and substitute into Eq. (5.3), which can then be written as 

V 2 h k0 + 2kh k0 = 2a 2 KSTl - U/ k V, (5.9) 

where the 'source' term on the r.h.s. is known. Considering, alternatively, Eq.(5.5), we get 

V 2 V + 3a(-n)'v = -\a 2 KV k ST k , (5.10) 

h k Q = -2a 2 nST k - \v k V. (5.11) 

B. Global gauge conditions and integral-constraint vectors for spherical universes 

Let us now consider the integral gauge conditions (4.57). We wish to elucidate their 
relation to Traschen's constraint vectors [37], [38]. An integral-constraint vector is defined 
by the relation 

/ ST^V^d^V = [ dY<iB l , (5.12) 

JV JdV 

in which V is (possibly a part of) a spacelike hypersurface, n a its normal, dV its 2- 
dimensional boundary; B l depends on h^ v and its derivatives and it vanishes if these are zero 



V 2 h k n - 2 



k - 2a 



In) 

a / 
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on dV; V n is gauge independent. Since does not depend on 5T^, Eq. (5.12) represents 
simple constraints on source perturbations. 

There exist 10 integral-constraint vectors in each of the FRW universes but 6 of them 
are just spatial Killing vectors. The other 4 are more interesting - Traschen and others 
considered their implications for microwave background anisotropies (see, e.g., [37], [39]). 
The time components of the 4 Traschen vectors are proportional to the scalars Q B [Eqs. 
(4.26), (4.27)], the spatial parts - to V l Q B . In a closed spherical universe 

V^ B) = {Q B ,a- l ixV i Q B ). (5.13) 

Applying (5.12) to the whole closed universe, it takes the form 



Q B ST o + ^qb St o 



I S 3 

Integrating by parts in the second term, we obtain 



53 L 



qb St o _ 1qb^ 5T o 



d {3) V = 0. (5.14) 



d (3) \/ = 0. (5.15) 



a 

In order to deduce the simplest constraints on the matter perturbations, Traschen et al. [37] , 
[38], [39] consider the synchronous gauge and, in addition, restrict "physics" in assuming 
vanishing pressure so that the synchronous coordinates can be chosen to be comoving with 
the fluid (in fact dust) since the flow is irrotational. Then the fluid velocity 14 = 9k\V l = 
gkiSU 1 = and 5Tf — (p + p)(hi + Vi) — in the synchronous gauge. Since the second 
integral in (5.14) vanishes in this case, the constraints imply just [see, e.g., (14) in 



/ Q B 5pd^V = 0. (5.16) 

In the Mach 1 gauge the constraints (5.14) and (5.15) have clear, simple consequences 
without any necessity to restrict physics. Since both the constraint equations (5.2) and (5.8) 
have on the l.h.s. the operator V 2 + 3 which has eigenfunctions Q B with zero eigenvalues 
(see Appendix C), it is evident that the sources on the r.h.s. must be orthogonal to the 
4-dimensional function space spanned by 4 harmonics Q B . Therefore, in closed universes 
the perturbations 5T{? of any type of matter have to satisfy separately the constraints 

/ Q B 5T°d ( Vv = 0, (5.17) 



and 



/ Q B V^V 3) V = 0, (5.18) 

JS3 
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the last being equivalent to 

I V i Q B 5Tfd {z) V = 0. (5.19) 
Js 3 

The same is true for 5T® = a _1 5T^° and 5T® = 5T®. Hence, Traschen's constraints (5.12), 
resp. (5.14), are indeed satisfied - in such a way that, in fact, both integrals in the constraints 
have to vanish separately. The constraints now become a straightforward consequence of 
the Einstein equations. This is not the case in the synchronous gauge where the constraint 
equations are coupled and there are more complicated equations for h^T^ = ^7ih l Tk and 
their derivatives, as can be seen from Eqs. (3.9) and (3.10) with hoo = hko = 0. 

The constraints (5.17) and (5.18) and the constraint field equations (5.2) and (5.8) also 
demonstrate, why our global gauge conditions (4.57) do not restrict physics. They just 
eliminate solutions of the homogeneous equations which, in any case, can be removed by 
gauge transformations generated by conformal Killing vectors. The gauge condition (5.1) 
implies that the same integral gauge constraint, satisfied for /i™ and V, is valid also for the 
spatial scalar hoo, as is also required in (4.57). As a consequence, from the field equation 
(5.4) another constraint, which has not been discussed by Traschen et ai, follows: 

/ Q B 5T"d (3 V = 0, (5.20) 
Js s 

the same for 5T™(= 5T™). Hence, in the Mach 1 gauge the whole picture of the Traschen-type 
constraints and our global integral gauge conditions is nicely symmetrical: all scalar pertur- 
bations in both the metric and the energy-momentum tensor, h^h^V = V fc/ig, ST™, <5T °, 
and V fe <5T°, are orthogonal to the 4-dimensional space spanned by harmonics Q B . 

C. Solutions of the field equations and local inertial frames 

We are interested in solutions for hoo,h i, and h™ when the sources are given in terms 
of 5T{f. These quantities determine local inertial frames. Hence, we wish to solve elliptic 
equations (5.2)-(5.4). Solutions can be given in terms of harmonics but these will not be 
considered here. However, several of these equations have been solved in literature in terms 
of Green's functions. The Green's functions for the equation 

V 2 $(^) + 3ife$(a: i ) = -2P(^), (5.21) 
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where x % = {x,9,ip} are the hyperspherical coordinates, are [39]: 
G s s(x,x') = -(1/4tt) 



cos 2ib ( 1 

— — — 1 — — cos ip 

sm ip \ 7i J 2tt 



G E3 (x,x') = -(l/4n)l, 



G H3 (x,x') = -(l/4n) 



cosh 2a 
sinh a 



— 2 cosh a 



(5.22) 
(5.23) 
(5.24) 



where cos ip = cos x cos x'+sin % sin cos 7 (k = +1),/ = l 2 +l' 2 — 211' cos 7 (k — 0), cosh a = 
cosh x cosh x' — sinh x sinh x' cos 7 (A; = —1), and cos 7 = cos 9 cos + sin 9 sin cos (</? — 
tp'). Here a is the geodesic distance under the metric between the 'source point' x H = 
{x', 9', ip'} and the 'field point' x l = {x, 9, <p}. The Green's functions satisfy the equations 



(V 2 + 3k)G(x,x r ) = [det(f i:j }—2S{x,x'), 



(5.25) 



where V 2 refers to the point x l , S(x',x') is the Dirac distribution. In terms of the Green's 
functions (5.22)-(5.24) the solution to Eq. (5.21) is given by 



$(ar) = ~ 2 J G{x,x')P(x')S 3) V. 



(5.26) 



Hence, given <5T ° and V fc 5T° we can determine h™ and V from Eqs. (5.2) and (5.8). h 00 can 
be determined from the gauge condition (5.1), or by solving (5.4) if 5T™ is known. Notice, 
however, that to determine /iqo from (5.1) we need to know also the time derivative h™. This 
can be found by taking the derivative of Eq. (5.2) and assuming that 5Tq is known. [5T® 
can be expressed from the Bianchi identity (3.25) in terms of 5T^ and h 00 .] The solutions 
for h™ can then be constructed. 

Knowing V from Eq. (5.8), we can determine hko from Eq. (5.9) or (5.11). In general, 
we need to find a Green's function bitensor G a b '(x, x') satisfying 



f lm V l V m G a b '(x,x') + 2kG a b '(x,x') = r^ 2 (x)5 b a 5(x,x' 



(5.27) 



in the case of Eq.(5.9) [analogously for Eq. (5.11)]. Then the solution for h k can be written 



as 



h ok = J G k b, (x,x')Sl A) (x')d^V'; (5.28) 
by the source S^ A \ A = 7, 77, the r.h.s. of Eq. (5.9), respectively (5.11), is denoted. 
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In a spatially flat universe, the easiest way is to write Eqs. (5.9) and (5.11) in Cartesian 
coordinates. Then (5.9) decouples into three Poisson equations for each h k0 , the Green's 
functions are standard, and the solutions are given by Poisson integrals over the source: 

hok(x\t) = / p-VLLtfx', SP = 2a 2 KbTl - -d k V. (5.29) 

J \X X \ o 

Equation (5.11) turns into three equations of the Yukawa-type, as noticed recently by Schmid 
[58]. Indeed, the l.h.s. of eq. (5.11) is of the form V 2 /iq — \ 2 (r])h^ , so two Green's functions 
are given by 

G(x, x) = — - : -, (5.30) 

47T \X — X'\ 

where 

\ 2 ( v ) = -4a(H/a)', \ 2 {t) = -4a 2 H. (5.31) 
Usually H < 0, so A is real. The well-behaved solution of Eq. (5.11) is thus 

h k = ~J ^'V)^^V, S™ = -2a 2 K5f k + \d k V. (5.32) 

For open universes, the properties of the Green's bitensor G a b ' solving Eq. (5.27) 
with k = — 1 have been studied by d'Eath [59]. In particular, it can be shown that such 
G a b ' exists which satisfies the boundary conditions at the source points and decays as 
exp[— 3d(x, x')\ as d(x,x') — > oo, with d(x,x') being the geodesic distance between the 
points x, x' under the metric of an open universe. In hyperspherical coordinates, d = a, 
where a is given below Eq. (5.24). In fact, it was d'Eath [59] who found the explicit form 
of the (scalar) Green's function (5.24), but the explicit form of the Green's bitensor for 
solving the equations for vector perturbations for k = — 1 does not seem to be known. The 
same is the case with spherical universes where only the Green's function (5.22) for the 
scalar equation (5.21) with k — +1 is known. Nevertheless, we can find explicit solutions 
for quite general classes of the vector perturbations also in case of k — ±1. 

1. Axisymmetric rotational perturbations 

Recently we solved Eqs. (5.9) and (5.11) for all odd-parity vector perturbations, i.e. those, 
for example, corresponding to rotational perturbations with axial symmetry [40], [41]. We 
decomposed perturbations in coordinates 9, ip on spheres only and assumed axial symme- 
try (spherical functions Yj m having m — 0). Since the backgrounds admit homogeneous, 
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isotropic foliations, nonsymmetric perturbations can be found from the axisymmetric ones 
[41]. Thus, we write in spherical coordinates of Eq. (1.2) 



oo 



i=i 

oo 



K = Y<[K^,r)\s\*eY lQ ,o, 



i=i 



(5.33) 



(5.34) 



where Yio,e = d$Yio, and 5Tjf may represent any perturbation. In the case of perfect fluid, 
the fluid angular velocity [cf. Eq. (3.16)] is Q = = dtp/dr/ = Q/a, and we write 



^ = -aV^Q,(t,r) sintflV 



(5.35) 



e=l 



Putting 

ho v = a 2 r 2 sin 2 9 uj i (t,r), (5.36) 

we have 

<5T°] i = a 2 (p + p)r 2 sin 2 9 (£, - (5.37) 

These perturbations are transverse: Vfe/i* = = VfcT fc , V — 0. Equations (5.9) and (5.11) 
become 

V 2 V + 2A;V = 2a 2 ^f °, (5.38) 

respectively, 

V 2 J# - 2[k - 2a(H/a)']h$ = -2a 2 K6T$. (5.39) 
The relation (3.15) now implies 

2 r 



8T° = -r 2 sin 2 9 



A ' _ 

0J 2 



(5.40) 



so Eq. (5.39) immediately follows from Eq. (5.38) and vice versa. Nevertheless, the equa- 
tions differ in the sense that in Eq. (5.38) 8T® is considered as a source, whereas in (5.39) 
the source is given by 5T£. <5T° determines (up to factor a 4 ) the density of the angular 
momentum - the perturbed Bianchi identities (3.23) imply the conservation law 



(a 4 6T°Y = 0. 



(5.41) 



On the other hand, STq determines the energy current. This is most apparent in the case 
of perfect fluid: ST® is given by Eqs. (5.34) and (5.37), while 



ST* = (p + p)V* = (p + p)Q = J2[STg}i(sm9)- 1 Y l0t g, 



(5.42) 
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where 



5T^(rj, r)j = (p + p)Q t . Substituting the expansions into Eq. (5.38), and using the 



orthogonality of sin^ F/ ,6> for different /'s, we obtain the "radial" equation for each I: 



Q 

Vl — kr 2 — (r 2 uji j 



dr 



+ r 2 Ul ~ Akui = 2a 2 «(p+P)(fii-^) = X^rti-Ui). 

(5.43) 

For I = 1 the perturbations correspond to the 'rigidly rotating spherical shells' in the FRW 
universes [12], [40]. Each sphere rotates with no shear but it expands/contracts with the 
background so that its angular velocity changes. For / > 2 the motion of the fluid is 
"toroidal" [41]. In the case of closed universes the Legendre equation which follows from Eq. 
(5.43) requires a special treatment. For example, for k = +1, functions ui(t,r) determining 
hop by (5.36) and (5.33) turn out to be (r = sinx) 

2, = 2A;(sin X )- 3 / 2 \v\ f - % - (ST^dx' + Q l 2 f - % - (S^dA , (5.44) 
[ Jo Wism i/Z x v J x Wism 1/Z x v J 

where Wi is the Wronskian of the functions V l 2 (x)i Q l 2 (x) which are derived from the deriva- 
tives of the appropriate Legendre functions with respect to their degree [41]. 

The properties of the solutions of Eq. (5.43) differ significantly according to whether we 
consider the right-hand side of (5.43), i.e., the angular momentum density <5T° as the source 
of uji, or we solve (5.43) for ui with fli given, i.e., with the angular velocity as the source. 

The rotation of inertial frames (2.30) is given by the angular velocity 



-u 1 



1 

2a 



oo oo 1 , 

.1=1 1=1 



(5.45) 



The complete solutions for uj\ for both 8T® and 8Tq given are determined in [41]. Roughly 
speaking, in flat and open universes the effects of torodial motions beyond the cosmological 
horizons are exponentially damped when the angular velocity of matter is given. For flat 
universes, this was first noticed by Schmid [58]. We shall see it occurs also for accelerations. 
However, these dragging effects are not damped when angular momenta are given as sources. 
In [40] we give the physical explanation. Since Eqs. (5.9) and (5.11) are elliptic equations, 
in both cases the inertial influences of "distant matter" are expressed instantaneously. 

We found toroidal perturbations to cause the rotation of local inertial frames by the 
angular velocity (5.44). Do they cause their acceleration? Since now only h 0(p ^ among 
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all h^, the only nonvanishing component of the acceleration (2.20) is 

1 



a* = 



W- (5-46) 



r 2 sin 2 9 

Substituting for h 0ip from Eqs. (5.33) and (5.36), and expressing the acceleration in the 
'background' frame, we find 

oo 

a (f) = ~ r $Z -( fl2 ^)'Ve. (5-47) 
i=i 

The acceleration vanishes at the axis of rotation. More interestingly, it vanishes everywhere 
in a static (Einstein) universe if the matter rotates uniformly. Indeed, the angular momen- 
tum density conservation law requires [a 5 (p + p)(u — fi)] = 0, which implies Co = and 
hence = for time-independent fl and a = constant. In the FRW universes the accel- 
eration (5.47) is nonvanishing, thus resembling the acceleration of the local inertial frames 
with respect to the static frames inside a collapsing, slowly rotating shell where particles 
at rest with respect to infinity experience the Euler acceleration, although the spacetime 
inside the shell is flat [32]. For all types of the FRW universes the solutions are of the form 
u>i = gi(r)a~ 3 (t), where gi(r) are explicitly given in terms of the integrals of the special 
functions mentioned above and the sources [5T<9]j. Hence, the accelerations (5.46) are of the 
form 



oo 

n, , = -(H/a 2 )r^2 gi (r)Y m , e . (5.48) 



i=i 



As an illustration, for k = and I = 1 perturbation we get 

a M = 2(H/a*)rsin9(^± + f ^dr'), (5.49) 

where J(< r) is the angular momentum inside r. With angular velocity Q considered as 
a source, u shows the exponential decline near the origin when the source is beyond the 
horizon [41]. As a consequence of Eq. (5.47) the acceleration behaves similarly. 



2. Perturbations of potential type 

In the example of toroidal perturbations, we had 5T® = V^V^ = Vfc<5T fc = V^/iq = h 00 = 0. 
As a second example, consider briefly the case in which these quantities may be nonvanishing 
but 14 and /i fe0 have a vanishing transverse part so that h^o = Vfc/i for some scalar h, and 
similarly for = V^w. Physically, such perturbations describe a change in the matter 
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density and a curl-free velocity field. No rotation of local inertial frames arises for such 
perturbations - the vorticity vector (2.29) vanishes for these "scalar perturbations". 

In order to determine the acceleration, we can use the gauge condition (5.1). It enables 
us to find /loo m terms of h™ and V. Both h™ and h™ can be determined from (5.2) and 
its time derivative. With STq given, the Green's functions (5.22)- (5.24) yield h™. In the 
simplest case of a flat universe, 

h n n = K [ (f SJ °\ d 3 x'. (5.50) 
J \x-x'\ 

The scalar V can be obtained by solving either Eq. (5.8) or Eq. (5.10). As with toroidal 
perturbations, when the angular momentum <5T° is prescribed, V = Vfc/ig will not be sup- 
pressed at the origin if 5T® occurs beyond a horizon. A suppression takes place if the velocity, 
or the energy current STq, is prescribed, as it corresponds to solving Eq. (5.10). In the flat 
universe, for example, Eq. (5.10) reads 

V 2 V + 3a 2 (£) V = ~a 2 KVi5fl (5.51) 

which is a Yukawa- type equation with the solution [A 2 (t) = —3a 2 (a/a)'] 

^n 2 K r ~ p -X\x-x'\ 

V^^f^S^f—. (5.52) 

If we start from Eq. (5.8) with 5T k given as a source, the solutions can be written in 
terms of the Green's functions (5.22) - (5.24). Taking d/dt of (5.8) and assuming 5T° k given, 
the same Green's functions will yield V . Since h ok = Vfc/i, V = V 2 /i, V = V 2 /i, we can find 
hko and hko by solving Laplace equations for h and h with V and V given. The solutions 
are unique up to an additive function of time which does not contribute to h ok . Knowing 
h k0 and h 00 = Icr 1 ^ — \ah^) from the gauge condition, the acceleration (2.20) of the local 
inertial frames can be determined in terms of the sources. 

Alternatively, we can start from Eq. (5.4) to determine directly h 00 in terms of the matter 
perturbations <5T ° and 5T™. Equation (5.4) for k = becomes just a Yukawa-type equation 
exactly in the form of Eq. (5.52). hoo exhibits an exponential suppression near the origin 
if the source 5T$ , ST™ occurs beyond the cosmological horizon. The suppression enters the 
formula (2.20) for the acceleration. However, the term /i fc will not be suppressed if the 
angular momentum 5T° k is prescribed. If the energy current 5Tq is considered as a source 
of hok, the total acceleration will be exponentially suppressed. In both cases, however, the 
acceleration is determined instantaneously. 
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D. The determination of local inertial frames 

We described how the accelerations and rotations of the local inertial frames can be 
determined explicitly in these specific examples in order to illustrate the general framework. 
Finally, let us discuss, within the Mach 1 gauge, a question of the uniqueness of the solutions 
of the field equations for general perturbations and of the resulting expressions for the 
acceleration and rotation of local inertial frames. The homogeneous equation corresponding 
to Eq. (5.3) for hko is identical to Eq.(4.23), the well-behaved solutions of which were 
analyzed in detail between Eqs. (4.23) and (4.33). They do not exist in H 3 ; in E 3 they 
describe the time-dependent linear combination of translations, which can be eliminated 
by requiring /i fe0 to decay at infinity. In S 3 , they correspond to the time-dependent linear 
combinations of 10 conformal Killing vectors. However, by imposing our integral gauge 
conditions (4.57), we dispose of the four conformal Killing vectors which are not Killing. 

Therefore, the complete general solution of Eq. (5.3) has for k — 1 the form 

6 

h 0k = h ok {inh) + ^ fA(t)ti A \ (5.53) 

A=l 

where h k (i n h) is a solution of the inhomogeneous equation (5.3) and /a(£) are arbitrary 
functions of time; £^ are 6 Killing vectors of S 3 describing rotations and quasitranslations 
(see Appendix C). Owing to our integral gauge conditions, Eq. (5.2) admits a unique 
solution; V = V fc /iofc is not affected by the Killing vectors in (5.53), and hence, the Mach 1 
gauge condition (5.1) determines a unique hoo- If we start from Eq. (5.4) to determine hoo 
directly in terms of 5T® and <5T™, we also arrive at a unique solution because the homogeneous 
equation corresponding to Eq. (5.4) coincides precisely with Eq. (4.38). Therefore, with <5T °, 
5T™, and angular momenta 5T® given, the accelerations and rotations of local inertial frames 
in closed universes are determined by formulas (2.20) and (2.31) only up to the freedom 
exhibited in Eq. (5.53). This freedom corresponds precisely to changing the coordinate 
system by the infinitesimal transformation in which £° = 0, C = Xm=i -^a(£)£^ 1 - Then 
A/ioo = A/i™ = Ahki = because ( l is a linear combination of the Killing vectors. The 
six spacelike Killing vectors generate motions which preserve the symmetries of the space. 
However, 

6 

A/^ = -aCk = -a ^mf, (5.54) 

A=l 



52 



which is equal to the additional term in Eq. (5.53) provided that /a(0 = —a(t)FA(t). The 
transformations 

6 

x * = x i + F A(t)Z {A)i (5.55) 

A=l 

with arbitrary coefficients lead to mutually accelerated frames. The accelerations have 

special forms when regarded as functions in space. Putting x l = constant in (5.55) we get 
cPx' l /dt 2 = J2a=i ^A{t)^ A ^ where are specific functions of x l (cf. Appendix C). An 
example of accelerations and rotations generated by this type of transformations with one 
translational and one rotational Killing vector is given in Eqs. (4.53) and (4.54) above. 

Consider now for simplicity just one rotational Killing vector and again the case of toroidal 
perturbations. The arbitrariness exhibited by the additional terms in Eq. (5.53) can then be 
seen distinctly. The homogeneous equation corresponding to Eq. (5.43) for I = 1 with the 
angular momentum given as a source is solved by u3( ) = cu^(t), where cD( ) is an arbitrary 
function of t. This implies [cf. Eq. 5.36] 

h 0tp = £( )a 2 sin 2 x sin 2 9 = UJ( )^\ (5.56) 

where = (0, 0, 0, 1) is the rotational Killing vector, which is a special case of Eq. (5.54). 
The transformation if' — ip — J U(p)(t)dt [a special case of (5.55)] would make the term (5.56) 
vanish. Since, however, in closed universes such an arbitrary 'integration constant' uJ(o){t) 
cannot be eliminated by boundary conditions, all frames with different a>(o) (t) are admitted. 
In this sense, only relative rotations of the local inertial frames can be determined if the 
angular momentum is considered as a source of their dragging. 

In the case of general perturbations of the spherical universes with the distributions of 
energy 5T$, angular momenta 5T® given, the freedom described by in Eq. (5.55) 

which preserves the symmetries of the space cannot be eliminated by boundary conditions. 
In this sense, only relative rotations and accelerations of the local inertial frames can be 
determined. The perturbations h 00 and h 0i , which imply these rotations and accelerations, 
are determined by the field equations and the Mach 1 gauge conditions instantaneously 
from appropriate averages over the distributions of <5Tq and 5Tf?. An explicit example is the 
expression (5.26) with the source term P ~ 5T® and the Green's function given in S 3 by Eq. 
(5.22) determining ft,™ and thus h 00 . Other example are the functions ui(t,r) in Eq. (5.44). 

This instantaneous determination of the local inertial frames by such averages, up to 
global rotations and accelerations given by the symmetries of the space, is the crucial feature 
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exhibiting the validity of Mach's principle in relativistic cosmology, at least in the first- 
order perturbation theory. The ability to describe the same physical situation using these 
differently rotating and accelerating frames is a consequence of the dynamics having a higher 
degree of symmetry than the realisation of the world in terms of the positions of actual 
bodies. The frame in which we choose to describe the motions is not of importance; what 
matters is the relative motions of the bodies, not that of the frame relative to the bodies. 

Finally, consider now energy currents 8Tq together with <5T ° as the sources. In case of the 
perfect-fluid perturbations we thus take the fluid velocity V k and Sp as the sources. When 
velocities and accelerations of "heavenly bodies" are given, the rotations and accelerations of 
the local inertial frames are determined uniquely in spherical universes. As we have shown 
below Eq. (4.38), the homogeneous equation corresponding to Eq. (5.10) for V = V^/ip 
admits only the trivial solution V — 0. Hence, the inhomogeneous Eq. (5.10) determines V 
uniquely when Vfc<5T fe is given. The same is true for (5.4) for h 00 because the homogeneous 
equation is the same. To determine perturbation ft o we need to know both Sp and Sp (resp. 
ST™). The gauge condition (5.1) then implies h™. Alternatively, we can solve (5.2) for h™ 
and ft™ by giving Sp and Sp, and extract ftoo from the gauge condition. Solutions for h™, h™ 
are unique due to our integral gauge conditions, so a unique h 00 can also be found. Finally, 
the homogeneous part of (5.11) for k — 1 with sources STq and V given reads 

h\ = 0. (5.57) 

This admits only = 0: multiply by hko and integrate by parts over closed space, 

- / rf^VihkoVjhiD^V = 2A(t) [ f M h ko h lo d^V, (5.58) 

J D JD 

where A(t) = 1 — 2a 2 (a/a)' = a 2 [|«(p + p) — H\. Since the integrands on both sides are 
spatial scalars, we can calculate them at each point by using f tj = diag (1, 1, 1). In this way 
we find that both are non-negative so that the only way how to satisfy (5.58) for A(t) > 
is by putting hko = 0. For standard models H < and indeed A(t) > 0. In the example 
of toroidal perturbations, the uniqueness of the solutions is reflected by the fact that Eq. 
(5.43) has unique solutions for given angular velocity Qi of matter. A purely time-dependent 
ui(p)i(t) does not solve (5.43) when Qi is prescribed and not the whole r.h.s. Qi — 

We thus arrive at another important aspect of Mach 's principle in relativistic cosmology: 
if the velocities, density and pressure perturbations of cosmic fluid are given, the (linearized) 
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field equations in a closed universe provide a unique determination of the rotations and 
accelerations of the local inertial frames. 

VI. CONCLUDING REMARKS 

Although this paper includes also items which have a review character, primarily it con- 
tains new developments: the analysis of accelerations and rotations of local inertial frames 
and of gyroscopes in perturbed FRW universes; the general forms of the perturbed Einstein 
field equations and Bianchi identities are formulated without gauge conditions, harmonics, or 
splittings, the motivation for and the analysis of the instantaneous, Machian-based gauges, 
including the integral gauge conditions and their relation to Traschen's integral constraints; 
and the manifestation of Mach's ideas in the framework of general linear cosmological pertur- 
bations of FRW universes. In particular, those who wish to study cosmological perturbation 
problems in position space, as advocated recently in Ref. [36], may find here useful relations 
not given before. Various specific perturbation problems can be attacked by applying the 
results presented here. We already used the formalism to investigate rotational and toroidal 
vector perturbations of FRW universes [40], [41], as mentioned and applied in Section V; 
there we also discussed vector perturbations of potential type. 

For given distributions of energy-momentum and angular momentum of matter sources, 
the rotations and accelerations of local inertial frames are uniquely given in the Machain 
gauges in open universes under suitable boundary conditions, whereas in closed universes 
they are determined up to motions generated by the Killing vectors, i.e., by symmetries of 
the background. They are determined uniquely also in closed universes if velocities, density, 
and pressure perturbations of cosmic fluid are given. As a consequence of the constraint 
equations and the choice of gauges which imply suitable slicing of perturbed universes, 
these inertial properties are determined instantaneously. In this sense Mach's principle is 
embodied in the cosmological linear perturbation theory. 

The dragging of inertial frames is an essentially global effect which, at least in linear 
perturbation theory, has to be seen as an instantaneous phenomenon. This was first demon- 
strated by Lindblom and Brill [60], who investigated rotational dragging by a slowly rotat- 
ing, massive spherical shell freely falling under its own gravity. We reconsidered the problem 
and explored its electromagnetic analogue [32]. The need to introduce a suitable coordinate 
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frame (the 'gauge') to describe the dragging is well illustrated inside the shell. Spacetime is 
flat there; no local geometrical (gauge invariant) perturbations occur. The time- dependent 
rotation of inertial frames is exhibited by considering the congruence of static observers, 
i.e. those who are at rest with respect to static observers at infinity. They play the role 
analogous to that of the cosmological observers in the present paper. They experience ac- 
celeration, and the congruence of their worldlines twists. Both quantities, characterizing 
their congruence, can be expressed in a covariant manner as in formulae (2.7) and (2.8) for 
cosmological observers. Also, massive, slowly rotating shells immersed in FRW universes 
were analyzed [61], [62], including their observational consequences on the appearance of 
sources behind the shells [62], [11]. In [33] we considered strong rectilinear dragging using 
exact conform static solutions of the Einstein-Maxwell equations with charged dust. 

A thorough nonlinear study of Mach's ideas within the framework of general relativity 
lies in the future. Quoting from the same source by which we started (see [1], p. 546), 
"Much must still be done to spell out the physics behind these equations [the initial-value 
equations] and to see this physics in action" . 
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APPENDIX A: PERTURBED FIELD EQUATIONS WITH COSMIC TIME t 

We write the perturbed FRW metric in the form (1.1) 

ds 2 = (g^ + h^) dx^dx v = dt 2 - a 2 (t)f H dx k dx l + h^dafdx" . (Al) 
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The background Christoffel symbols are 

rS = HST, T° kl = -Hg kl , fj5 = f mn (%/0n - \d n f k ^j ; (A2) 

hereafter the symmetrization brackets ( ) include the factor |, as do antisymmetrization ones 
[ ]; H = a I a is the Hubble "constant". The non- vanishing components of the background 
Einstein equations, = B£ - \8^R = kT» + AS^, read 

G° = 3^ + ^= K p + A, G[ = # k (± + 3H 2 + 2H^=-( K p-A)$ k , (A3) 

with k = 87rG/c 4 , k — 0, ±1 denoting the curvature index, and A cosmological constant; 
the background energy-momentum tensor T" of perfect fluid is given by (2.32). The indices 
of are raised or lowered with g^ and g^; thus h® = h 00 , = g kl h i = — ^f kl h i, etc. 
No spatial index is ever displaced with f ki alone. We introduce the covariant derivative, 

v fc zc = d k h™ + r$X- (A4) 

The background curvature tensor of spatial sections t = constant is VJ ks i = k (5 r s f k i — ^Ifks) 

and the Ricci 3-tensor TZ k i = 2kf k i. Useful identities are (Vfc, = VfcV,, V fc = f kl Vi, 
V 2 = / H V fc Vz): 

2k 2k 

(V fc , - V lk )V l = -2kf kl V l = —g kl V l = —V k , (A5) 

a 1 a 1 

Vfe VV fc = V 2 {V k V k ) + 2k(V k V k ) , (A6) 

(Vfc,- V lk )h l = 2^h ko , (A7) 

(V fei - V lk )h n m = h r m n n rkl - h n r K r mkl = 2k (Sf k h\ } - hfrSfi f rm . (A8) 
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The perturbed Einstein equations, 5G^ = k,5T^, are expressed in terms of h° k , and h l k . 
In this 'mixed' form A does not appear. The left-hand sides, 8G^ , read as follows: 



SG° = -\v rs {g rn K - f s K) - ^h n n - 2H {^-Hhl - l -K + v n K 

SG° k = \ Vt (h{ - 8{hl) + g rs (V kr h° s - V r(k h° s) ) + HVX , 



SG k = g kl 



6G°-2 i^-H )h? 



(A9) 
(A10) 
(AH) 



+ \ v« {g ln K - g lm K) - \g rs v rs (K - «) + \ v rs (g mr K - s?h r 
+ g ml (v (Jfc fc?) + HV {k hty - 5rg rs (v r h° a + HV r h°) 

- \ {g ml Vi k h° - 5^g rs V rs h° ) - 8? [Hh° + (2H + 3tf 2 ) h° . 



(A12) 



If 8T£ is a perfect-fluid perturbation then the right-hand side is given by [see Eqs. (2.35)] 

2 



5T ° = 8p , 8T° k = (p + p)(V k + hl) = -[- 2 -H)(V k + hi) , 

k \ci j 



8T k = (p + p)V k 



2,k 



k a" 



-H)V k , 8T{ = -8 l k 8p. 



In the last equations we used the relation 



k \a 



(A13) 



(A14) 



which follows from the background Einstein equations (A3) for all p, p, k, A. 
The relations between various h^s and h^s used in the main text are 



h 



oo 



hnn , ht 



01 



a hi 



01 



hki — a 2 h k i 



ho — hoo 



h l = —a 1 h l , h® = ah 



oi , h[ — -h[ , 
a- 2 h kl . 



(A15) 



h 00 = h 00 , h 01 = -a- 1 ^ , h kl 
Eqs. (A4) and (A6) hold also for h™ and V k , but Eqs. (A5) and (A7) take the form 

(Vh - V lk )V l = -2kf kl V l = -2kV k , (A16) 
(V fc , - V lk )h l = -2kf kl h l = -2kh ok , (A17) 

and Eq. (A8) becomes 

(Vh - Vik)h n m = 2k (5? k h l]m - hlf l]m ) . (A18) 
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APPENDIX B: GAUGE TRANSFORMATIONS OF PERTURBATIONS 



As a consequence of infinitesimal transformations 



x°^x /0 = x° + C(x), 



x i -> x H = x i + C(x), 



(Bl) 



:b2) 



we find the following changes of various metric components under the change of gauge (notice 
that AQ = Q - Q' for any Q): 

Ah 00 = 2(° = Ah 00 , (B3) 
Ah ol = d t (° - a 2 Ci = aAh oh (0 = M k ) 



Ah 



ki 



-2a 2 



Ah" = 2a 2 



V(fcO) + -fklC 

a 



V n C + 3-C C 



a 2 Ah kt , 



-Ah" 



(B4) 
(B5) 

(B6) 



Similarly, the perturbations of the energy-momentum tensor components change under the 
transformations (Bl) and (B2) as follows: 



A5T ° 
A5Tt 
A 5T* 
AST l k 



pC° = A5f °, 
(p + p)d k ( = aA5T° k , 
-(p + PK^a-'AST*, 
-p(% = A6Tl, 



(B7) 
(B8) 
(B9) 
(BIO) 



where we substituted from Eq. (2.32) for the background values of T^. In particular, in the 
fluid case 



A6p = p(°, A5p = p(°, 
A5U = -(° = -A5U , A5U m = -C, A5U m = (° m , 
AV m = -C = a- l AV m , AV m = a 2 ( m = -aAV m . 



(Bll) 
(B12) 
(B13) 



Let us emphasize that the above results for the changes of both h^s and h^s are expressed 
in x M = {t,x l ) coordinates. In x^ = {rj,x l ) coordinates we find, for example, 



^oo ~~ ^-oo — A/ioo — ^-oo ~~ 2^° — /iqo ~~ 2 



^p-a- 1 = h 00 - 2HC° - 2^f. (B14) 
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APPENDIX C: KILLING AND CONFORMAL KILLING VECTORS ON THE 
FRW BACKGROUNDS 



Killing vectors are also conformal Killing vectors but here we call "conformal Killing 
vectors" - sometimes more explicitly "proper" conformal Killing vectors - those which are 
not Killing vectors. All these vectors on the FRW backgrounds are well-known. Since, 
however, we did not find listed all of them in a transparent manner in one place, we give 
them here. Their relation to the scalar and vector harmonics will also be elucidated. The 
3-dimensional spatial vectors are frequently used in the main text. There exists an extensive 
literature on the harmonics in S* 3 and if 3 , see., e.g., [57] and [63]. 

1. Killing and conformal Killing 3- vectors in E z (k = 0), S 3 (k = 1) and H 3 (k = — 1) 

The standard Killing equation 

V fc & + Vi& = (CI) 

in the FRW 3-backgrounds with curvature tensor (A8) can be written in an equivalent form 

V 2 & + 2A& = 0. (C2) 

The Killing vectors have their simplest form in the coordinates x m in which the metric 
is (see e.g. [3], Ch.13), 



ds 2 = dt 2 - a 2 it) 
where 



hi + 



dx k dx\ (C3) 



r 2 = (x 1 ) 2 + (x 2 ) 2 + (:r 3 ) 2 . (C4) 
The 3 (quasi)translational Killing vectors are given by 

Z w = VT=kr*6*,, J =1,2, 3, (C5) 

and the 3 rotational Killing vectors by 

£(./)< = e iJ x i^ ( C6 ) 

Here e is the usual permutation symbol, £i 23 = +1, with indices moved by 5^, resp. 5 lk . 
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The 4 conformal Killing vectors are also simple in the x l coordinates: 

ft = Vl - IferV, (C7) 
^(J)i = 5 iJ _ kx i x J^ k = ±1> ( C8 ) 

= _ x i x J^ k = Q ( CQ ) 

The same Killing and conformal Killing vectors in hyperspherical coordinates are more 
complicated but are directly connected with better known forms of (hyper) spherical har- 
monics. Here we denote r = sinx (k — +1), r = x {k — 0), r = sinhx (k = —1) and, 
correspondingly, r' = cosx, 1, coshx- As in the main text, we denote the six Killing vectors 
by ^ A ^,A = 1,...,6, and the four conformal Killing vectors by ip^ l ,A = 1,...,4. The 3 
(quasi)translational Killing vectors read 

£(i)» _ ( s j n cos ^ f'f- 1 C os 9 cos ip, —r'r^ 1 sin ip/ sin 9), 

£( 2 ) 1 = (sin 9 sin (p^r'r^ 1 cos 9 sin ip, r'r^ 1 cos ip / sin 0) , (CIO) 

f( 3 )< = (cos#, -r'r -1 sin 0, 0). 

The 3 rotational Killing vectors (C6) turn into the same forms independent of k: 

£W = (0,- sin <p,- cot 9 cos <p), £ (5)i = (0,cos^, -cotflsin^), f (6)i = (0, 0, 1). (Cll) 
The 'dilatation' conformal Killing vector (C7) is for all k simply given by 

^ = (r,0,0), (C12) 

whereas the other three conformal Killing vectors read for k — ±1, 

^(2)1 _ (/sin^cosy^r^ 1 cos # cosy?, — r -1 simp/ sin 9), 

^,(3)j _ ^ r > g j n ^ g j n ^ r -i cos ^ sin y?, r _1 cos </?/ sin#), (C13) 

^,(4)i = (r'sin^^-r-igm^ o), 

and for = 0, 

^,(2)1 _ Ir 2 (— sin cos </?, r _1 cos0cosy9, — r _1 siny^/ sin0), 

,0(3> _ l r 2 (_ s i n # s i n <^ r - x C os0 sin y2,r _1 cosy?/ sin6>), (C14) 

^,(4)i = i r 2^_ cog ^ _ r -i sin ^ 
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2. Scalar harmonics in S 3 and H 3 



In S 3 , the scalar harmonics Qum with L > I > 0, L,l integers, m = — /, +/, satisfy 

V 2 Q L i m + L(L + 2)Q Llm = 0. (C15) 
In the normalized form they read 

Qum = yfN^-^=P^\cos X )Yi m {0,ip), (C16) 

where Nli = <X+1 ( ^^ +1 - >! , are Legendre functions of the first kind, and Y tm are the usual 
spherical harmonics. In H 3 , the harmonics are Q\i m , A > is continuous and must be real 
for square integrability; they satisfy 

V 2 Qxim + (A 2 + l)Qxim = 0. (C17) 

The normalized form is 

Qxi m = ^i^=P~f^ (cosh X )Y lm (0, if), (C18) 
Vsmhx 5 

AT Ai = A 2 (A 2 + l)(A 2 + 2 2 )...(A 2 + / 2 ). 

It can be easily seen that, for L — 1, I — 0, m — and L — 1, I — 1, m — —1, 0, +1, the 
expression (C16) leads, up to multiplicative constants, to the four functions given in Eq. 
(4.26), whereas nonintegrable harmonics (C18) for A = 2i, I = 0, m = and A = 2i, I — 1, 
m = —1,0, +1 imply Eq. (4.27). Their gradients (4.25) yield the covariant components of 
the conformal Killing vectors (C12) and (C13). 

Let us remark that for k = ±1 the translational Killing vectors are not gradients of scalars 
(as they are for k — 0). They are proportional to the vector spherical harmonics with even 
parity and L — 1, I — 1 (for k = ±1) and A = 2i, I — 1 (for k = —1). The rotational Killing 
vectors are proportional to vector harmonics with odd parity and L = 1, / = 1, respectively, 
A = 2i, I = 1. 
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APPENDIX D: FIELD EQUATIONS AND SOLUTIONS IN OTHER GAUGES 
1. Mach 2 gauge 

The gauge conditions T k = Vih l Tk = and V 2 h™ + 3kh n n = [cf. (4.3) and (4.9)] simplify 
the field equation (3.10) into the relation 

-2WC = a 2 K5T°, (Dl) 

from which /C can be expressed and V 2 /C needed in the following step can easily be 
calculated. Applying V fc to Eq. (3.11) we obtain the equation for V = Vih l , 

V 2 P + 3kV = la 2 K(V k 5T° + -i-V 2 <5T °). (D2) 
2 3rt 

Solving for V and substituting back into Eq. (3.10) we get the elliptic equation for h ok : 

V% k + 2kh 0k = 2a 2 K{bfl + ^V fe 5f °) - X - V k V. (D3) 

From Eq. (3.12) h 00 can be determined. /C can be expressed in terms of 5T® from relation 
(Dl), and the last the term on the r.h.s. involving the time derivative (a 2 7i _1 5TQ )' can be 
calculated by employing the perturbed Bianchi identities, Eq. (3.22). Eq. (3.12) becomes 

V% + 3kh 00 = ^(V k 5T° k - ^f °). (D4) 
ri ri 

From the elliptic equations (D2)-(D4) the metric perturbations h m , h ok follow instanta- 
neously if the sources 5T$ and (resp. STq) are given. 



2. Mach 3 gauge 

Together with T k = Vih l Tk = it is now assumed V 2 V + 3kV = 0, or simply V = V;/i = 
0. Applying V fc to Eq. (3.11), one gets 

V 2 /C = ^a 2 KV k 5Tl (D5) 

With /C known we obtain the elliptic equation for from Eq. (3.10): 

V 2 h n n + 3kh n n = 3a 2 K5T% + 6WC. (D6) 
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Next, we make the time-derivative of (D5) and substitute for V k 5T' k from the perturbed 
Bianchi identities. Applying then V 2 on Eq. (3.12) we arrive at the elliptic equation for h 00 : 

V 2 (V 2 /^o + Zkhao) = a 2 K [v 2 (<5T ° - ST%) + 3HV k 5T° + 3V fe V m 5T fc m ] . (D7) 

There is another simple elliptic equation satisfied by the quantity x — ^oo — Taking 
V 2 of Eq. (D6) and regarding Eq. (D5), we combine it with Eq. (D7) to obtain 

V 2 (V 2 X + 3A; X ) = 3a 2 «(V fe V^T< - l -V 2 5f:). (D8) 

The quantity x appears directly also in Eq. (3.13) for the spatial components SG l k . Applying 
V fc Vz on this equation, one arrives again at Eq. (D8) above. 

With V known (V = in the simplest choice of the Mach 3* gauge) and K. determined 
from Eq. (D5), the constraint equation (3.11) becomes a simple elliptic equation for h k0 : 

V% + 2kh k0 = 2a\6T° - ^V k V - ^V fc /C. (D9) 

The equations for h 00: h ok in both Mach 2 and 3 gauge are elliptic. Their form is very 
similar to the equations in the Mach 1 gauge. We can solve them by the same methods. 

3. Generalized Lorenz-de Donder gauge 

We start from the gauge conditions (4.15), expressed explicitly in Eqs. (4.16) and (4.17). 
Now in general T k = Vih l Tk ^ 0, and also V = Vih l and K. are nonvanishing. Nevertheless, 
the field equations (3.10)-(3.13) can be rewritten into a quite telling form. Denoting h 00 = 
(p, = ip, V — Vih l , we arrive at the following system: 

VV -<p"- 2Ha 3 - 6a (^j ip - AHV = a 2 n(5f° - 57%), (D10) 

VV - V" + ^ - 2Ha 3 (J^J - 2a (^j up - ^HV = a 2 n(5f° + ±6T%), (Dll) 

V 2 h 0k - Kk + 2kh ok - 4(Hh 0k )' - HV k (ip + 3ip) = 2a 2 ^f fc °, (D12) 
V 2 p - V" + m> - 4(HV)' - HV 2 (<p + 3i/j) = 2a 2 reV fe <5T°, (D13) 

-V 2 ^+^^+2^ fe +2^^+4W(/^V (m ^ )0 -^4V^) = 2a 2 K{5f l k - l -5 l k 5f%). (D14) 
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All equations now have the character of hyperbolic generalized wave equations. The metric 
perturbations are not determined instantaneously in terms of the sources 5Tjf. Although the 
main parts of the equations are given by the standard wave operators V 2 — d 2 /d^ 2 , there 
are terms involving lower derivatives of the metric perturbations which make the system 
coupled. These equations may turn out to be useful in cosmology as the standard harmonic 
gauge is in the post-Minkowskian approximations to general relativity. 
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